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A  general  mathematical  investigation  of  the  instability  prior  to  physical  breakup 
of  the  annular  liquid  fuel  flowing  from  the  airblast  nozzle  using  airswirlers  in  gas  turbine 
combustors  is  presented.  This  investigation  seeks  to  give  a  clear  understanding  of  the 
growth  of  instabilities  in  an  annular  liquid  fuel  sheet.  This  understanding  will  help  guide 
atomizer  design  efforts  to  achieve  improved  spray  quality  and  emissions  reductions. 

The  annular  fuel  sheet  flow  under  consideration  is  assumed  to  be  inviscid  and 
incompressible.  The  mathematical  formulation  applied  linearized  theory  to  the  Euler 
equation.  The  mathematical  model  derived  is  a  single  second  order  ordinary  differential 
equation  using  modified  Bessel  functions  of  the  first  and  second  kind. 
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Experimental  data  necessary  for  the  mathematical  investigation  are  provided  by 
the  University  of  California,  Irvine,  Combustion  Laboratory.  The  mathematical  model  is 
solved  for  the  growth  rate  of  a  disturbance  of  helical  fuel  sheet  using  a  computer  code. 
Variable  parameters  are  axial  wave  number,  tangential  wave  number,  surface  tension, 
axial  velocity,  angular  velocity,  density  and  the  radius  at  the  interface  in  the  mathematical 
model.  The  analytical  solutions  of  axial,  rotating,  and  helical  flow  annular  jets  were 
derived.  When  the  growth  rate  for  an  axial  flow  annular  jet  was  investigated,  angular 
velocity  and  tangential  wave  number  were  not  considered.  The  computer  code  was 
written  in  Fortran,  linked  to  IMSL  routines  and  executed  on  a  CRAYC90. 

The  instability  growth  rate  showed  a  difference  between  rotating  flow  annular  jet 
and  helical  flow  annular  jet.  The  growth  rate  for  a  helical  flow  annular  sheet  is  much 
larger  at  small  axial  wave  numbers  than  that  of  the  axial  flow  annular  sheet.  Also,  the 
growth  rate  for  the  helical  sheet  increases  with  the  angular  wave  number  while  the  growth 
rate  for  the  helical  sheet  increases  only  at  very  small  axial  wave  numbers  and  then 
decreases.  These  results  seem  to  be  enough  to  explain  the  phenomenon  of  sudden  fuel 
sheet  breakup  just  after  the  nozzle  tip.  The  effects  of  thickness  decrease  and  angular 
velocity  increase  are  relatively  strong  in  comparison  with  other  effects.  They  are  the  best 
candidates  to  be  used  in  the  design  of  atomizers. 
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CHAPTER  1 
INTRODUCTION 

1.1  Physical  Considerations  and  Objectives 
When  liquid  fuel  is  injected  from  the  annular  nozzle  injector  into  a  gas  turbine 
burner  section,  the  fuel  swirls  circumferentially  as  it  proceeds  downstream,  in  an  annular 
sheet  with  air  surrounding  it.  Figure  1  shows  the  location  of  the  annular  nozzle  injector 
in  the  typical  gas  turbine  engine.  Figures  2  and  3  show  the  atomizer  head  cross  section 
and  atomizer  head  cutaway.  Figure  4  shows  an  actual  photograph  of  annular  fuel  sheet 
disintegration  just  after  the  nozzle  tip  in  the  experiment  by  McDonell  and  Samuelsen 
(1988). 

Aerodynamic  forces  on  the  inside  and  outside  of  the  fuel  sheet  produce 
instabilities  that  eventually  cause  breakup  of  the  fuel  sheet.  Also,  there  are  the 
instabilities  due  to  different  centrifugal  forces  of  the  liquid  and  the  air.  Since  both  the 
air  and  the  fuel  have  high  velocities  in  the  axial  direction  and  circumferential  direction, 
even  disturbances  of  very  short  wavelengths  exist  in  both  the  axial  direction  and  the 
circumferential  direction.  These  instabilities  are  known  as  Kelvin-Helmholtz 
Instabilities.  In  addition,  centrifugal  instabilities  occur  because  the  liquid  fuel  is  much 
denser  than  the  air  outside  the  fuel  sheet.  The  outer  air  being  much  less  dense  moves 
radially  inward  while  the  denser  fuel  moves  radially  outward. 
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Figure  1.  Typical  Aircraft  Gas  Turbine  Engine 


outer  airswirling  flow  ■ 


Figure  2.  Atomizer  Head  Cross  Section;  Exit  Plane 


^>  Outer  air  flow 


~J>   Outer  air  flow 
Figure  3.  Atomizer  Head  Cutaway 
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Figure  4.  Schematic  of  spray  structure  in  experiment  at  UCICL 


The  combination  of  Kelvin-Helmholtz  Instability  and  centrifugal  instability  is 
expected  to  have  an  impact  on  the  finely  atomized  liquid  fuel  just  beyond  the  nozzle  tip 
in  the  gas  turbine  combustor.  Most  importantly  the  centrifugal  instability  is  expected  to 
dominate  the  sudden  breakup  of  the  liquid  jet.  To  validate  the  strong  impact  of 
centrifugal  instability  on  the  finely  atomized  annular  liquid  fuel,  the  mathematical 
formulation  of  the  helical  disturbance  must  be  derived.  The  mathematical  expression 
used  to  solve  the  problem  of  those  instabilities  actually  consists  of  a  system  of  nonlinear 
partial  differential  equations.  The  usual  approach  to  an  instability  problem,  assumes  that 
for  small  disturbances  the  equations  may  be  linearized.  The  terms  quadratic  or  higher  in 
the  disturbances  and  their  derivatives  are  neglected.  The  dispersion  relationship  will  then 
be  derived  from  the  application  of  interfacial  conditions  and  boundary  conditions  to  the 
mathematically  linearized  equation. 

Objectives  of  this  research  are  1)  to  describe  the  mathematical  formulation  which 
describes  the  growth  of  the  instability  for  an  annular  sheet,  2)  to  solve  the  mathematical 
formulation  using  a  computer  code,  3)  to  validate  the  computer  code  by  comparison  to 
other  known  solutions,  and  4)  to  study  quantitatively  all  the  effects  of  fuel  surface 
tension,  density,  fuel  sheet  thickness  and  air  swirling  on  the  growth  of  instability. 

1.2  Literature  Review 


The  process  of  breaking  up  a  liquid  jet  into  droplets  of  diameters  much  smaller 
than  the  jet  diameter  is  called  jet  atomization.   This  process  has  various  applications, 


including  the  liquid  fuel  injection  in  internal  combustion  engines  and  gas  turbine  engines. 
A  good  knowledge  of  the  fundamental  mechanism  of  atomization  is  essential  for 
increasing  combustion  efficiency  and  reducing  the  formation  of  pollutants.  Despite  its 
practical  importance,  the  mechanism  of  atomization  is  not  yet  well  understood. 

Advanced  gas  turbine  engines  are  dependent  upon  the  development  of  advanced 
fuel  injector  systems  that  will  provide  an  efficient  and  very  uniform  combustion  system 
over  a  wide  operating  range.  Emissions  of  advanced  gas  turbine  engines  have  been  the 
subject  of  many  studies  and  programs.  Of  particular  importance  are  the  oxides  of 
nitrogen  (NOx)  which  are  produced  in  advanced-cycle  engines  that  operate  at  high 
pressures  and  temperatures. 

For  combustion  to  be  effective  and  environmentally  acceptable,  the  fuel  must  be 
atomized  well.  If  good  atomization  is  not  obtained,  nonuniform  burning  will  occur. 
Localized  hot  spots  at  near-stoichiometric  conditions  cause  the  highest  temperatures  and 
produce  the  greatest  amount  of  NOx.  Reducing  the  NOx  emissions  by  reducing  the 
reaction  temperature  could  adversely  affect  the  combustion  stability  and  could  increase 
emissions  of  carbon  monoxide(CO)  and  unburned  hydrocarbons,  especially  at  low-power 
operation.  Therefore  a  means  of  controlling  the  reaction  zone  stoichiometry  over  the 
entire  operating  range  is  required  to  minimize  NOx  production  while  maintaining  good 
combustor  efficiency.  Uniform  fuel  injection  and  mixing  are  essential  for  achieving  these 
conditions. 

Soot  formation  is  also  associated  with  poor  fuel  atomization.  Soot  is  not  an 
equilibrium  product  of  combustion,  and  therefore  its  formation  is  influenced  as  much  by 
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the  physical  processes  of  atomization,  evaporation,  and  fuel-air  mixing  as  by  reaction 
kinematics.  Soot  is  generally  produced  anywhere  that  fuel-air  mixing  is  inadequate. 
Problems  of  concern  associated  with  soot  formation  are  higher  liner  temperatures  due  to 
increased  radiative  heat  transfer,  impingement  of  carbon  on  metal  surfaces  resulting  in 
erosion  and  reduced  equipment  lifetimes,  and  distortion  of  fuel  spray  distribution  because 
of  carbon  deposits,  leading  to  hot  spots  and  increased  NOx  formation. 

The  study  for  the  atomization  of  liquid  jet  has  a  long  history.  Rayleigh(1879)  first 
studied  the  phenomenon  for  the  breakup  of  liquid  jets  into  drops  mathematically  by 
treating  the  jet  as  a  doubly  infinite  parallel  flow  and  determining  when  that  flow  becomes 
unstable  because  of  temporally  growing  instability  waves.  It  was  natural  for  the 
experimentalists  to  approach  this  problem  by  imposing  controlled  harmonic  (usually 
acoustic)  disturbances  of  a  single  frequency  on  their  flows  (either  from  within  or  from 
outside  their  nozzles)  which  then  generate  spatially  amplifying  waves  that  are  closely 
related  to  the  temporally  growing  instability  waves  of  Rayleigh.  Castleman  (1931)  did 
some  of  the  early  work  in  this  area.  He  postulated  that  jet  breakup  is  due  to  aerodynamic 
interaction  between  the  liquid  and  gas  phases  on  a  moving  jet.  He  argued  that  this 
interaction  leads  to  the  formation  of  ligaments  or  liquid  columns  which  have  one 
end  still  attached  to  the  surface  of  the  jet,  and,  once  formed,  these  ligaments  rupture 
according  to  Rayleigh's  classical  breakup  theory.  He  proposed  that  the  formation  of  the 
ligaments  was  due  to  the  relative  gas-liquid  motion.  He  observed  from  the  experimental 
results  that  the  size  of  the  ligaments,  and  hence  their  lifetimes,  decreases  producing 
smaller  drops  as  the  relative  velocity  is  increased. 
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Schweitzer  (1937)  noted  the  existence  of  a  radial  component  of  velocity  in 
turbulent  pipe  flow  and  argued  that  this  radial  flow  could  cause  the  disruption  of  the  jet 
surface  as  soon  as  the  jet  emerged  from  the  nozzle.  He  separated  the  effect  of  the 
surrounding  gas  in  a  series  of  experiments  by  injecting  liquid  jets  into  evacuated 
chambers.  He  found  that  jet  atomization  could  not  be  achieved  unless  the  jet  was  injected 
into  a  compressed  gas.  This  led  him  to  reason  that  the  effect  of  turbulence  in  the  liquid 
flow  is  to  disturb  the  jet-gas  interface  and  to  augment  the  aerodynamic  jet  breakup 
effects.  He  proposed  that  if  the  emerging  jet  were  truly  laminar,  no  immediate  breakup 
would  occur  and  hence  liquid  turbulence  is  essential  to  the  jet  breakup  mechanism. 

Holroyd  (1933)  used  dimensional  reasoning  to  postulate  that  a  balance  between 
centrifugal  forces  (arising  from  angular  turbulent  velocity  components)  and  the  opposing 
surface  tension  may  account  for  the  jet  breakup  process.  Liquid  supply  pressure 
oscillations  have  been  noted  by  Giffen  and  Muraszew  (1953)  to  have  an  effect  on  the 
outcome  of  the  jet  breakup  process. 

Ranz  (1956)  analyzed  the  forces  involved  and  the  energy  expended  in  the 
ligament  model  of  breakup.  He  considered  the  ligament  to  be  stretched  by  the  action  of 
the  relative  motion  of  the  gas  and  liquid  flows  and  found  that  the  surface  energy  creation 
rate  was  comparable  to  the  viscous  dissipation  rate  and  that  the  surface  tension  forces  and 
viscous  forces  opposing  the  stretching  process  were  approximately  equal.  He  postulated 
that  the  ligament  diameters  were  related  to  the  wavelengths  of  unstable  aerodynamically 
induced  surface  waves  growing  on  the  core  of  the  liquid  jet.  He  proposed  that  the 
estimates  of  Taylor  (1940),  which  were  derived  from  considerations  of  the  growth  of 


infinitesimal  surface  waves,  could  be  used  for  a  quantitative  analysis  of  the  jet  breakup 
process. 

Rupe  (1962)  observed  that  high  velocity  laminar  liquid  jets  could  be  more 
unstable  than  fully  developed  turbulent  jets.  This  was  contrary  to  previous  information  in 
the  literature.  He  believed  that  this  behavior  was  due  to  changes  in  the  fully  developed 
laminar  velocity  profile  within  the  liquid  jet  as  the  jet  penetrates  the  chamber  gas.  This 
velocity  profile  relaxation  process  occurs  by  a  mechanism  of  momentum  transfer  between 
adjacent  liquid  layers  once  the  constraint  of  the  nozzle  wall  is  removed  at  the  nozzle  exit. 
The  kinetic  energy  per  unit  mass  of  the  discharged  liquid  is  dependent  on  the  details  of 
the  initial  fluid  velocity  profile  at  the  nozzle  exit.  In  pipe  flow,  with  a  parabolic  velocity 
profile,  the  kinetic  energy  of  the  fluid  is  twice  what  it  would  be  if  the  fluid  had  a  uniform 
velocity  profile  and  had  the  same  average  velocity.  In  the  course  of  the  passage  of  the  jet 
into  the  chamber  gas,  the  profile  relaxes  to  a  nearly  flat  profile.  He  suggested  that  the 
excess  kinetic  energy  is  redistributed  and  converted  into  potential  energy  within  the  jet. 
He  described  that  this  produces  a  radial  pressure  gradient,  and  hence  radial  velocities 
within  the  jet  since  the  kinetic  energy  is  concentrated  along  the  core  of  the  jet.  He 
postulated  that  these  radial  velocity  components  could  disrupt  the  jet,  lead  to  the 
formation  of  ligaments  which  could  then  be  broken  up  by  aerodynamic  forces. 

Generation  of  ripples  by  wind  blowing  over  a  viscous  fluid  was  investigated  by 
Batchelor  (1963)  with  linear  stability  analysis. 

Dombrowski,  N.  and  Johns,  W.  R.  (1963)  derived  the  linearized  formulations  of 
the  two-dimensional  planar  sheet  nozzle.  They  considered  a  sheet  which  is  thinning  with 
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time  and  moving  in  the  x-direction  with  a  constant  velocity  through  the  stationary  gas. 
They  obtained  the  equation  of  motion  of  the  neutral  axis  midway  between  the  two  gas  / 
liquid  interfaces  by  interrelating  the  forces  caused  by  gas  pressure  surface  tension,  liquid 
inertia  and  dynamic  instability  and  disintegration  of  viscous  planar  sheets  produced  by 
spray  viscosity.  They  described  that  the  displacement  of  the  liquid  sheet  is  proportional 
to  the  exponential  growth  with  time.  In  their  analysis,  the  disturbance  of  wave  amplitude 
is  higher  in  inviscid  sheet  than  in  viscous  sheet.  They  described  that  the  growth  of 
disturbance  is  increased  and  then  decreased  as  wave  number  is  increased  and  the  growth 
rate  of  disturbance  is  increased  as  the  thickness  of  sheet  is  decreased.  They  derived  the 
relation  between  sheet  thickness  and  the  diameter  of  the  ligament  which  is  formed  when 
disintegration  occurs  through  aerodynamic  wave  motion.  They  described  that  waves 
grow  on  the  sheet  until  they  reach  a  critical  amplitude.  He  suggested  that  tears  occur  in 
the  crests  and  troughs,  and  fragments  of  sheet,  corresponding  to  one-half  wavelength,  are 
broken  off.  They  also  described  that  the  fragments  contract  by  surface  tension  into 
unstable  ligaments  which  subsequently  break  down  into  drops. 

Shkadov  (1970)  investigated  the  effect  of  changes  in  the  interface  tangential 
stresses  in  a  stability  analysis.  He  considered  the  planar  problem  and  described  the  liquid 
flow  by  the  Boundary  Layer  Equations,  with  a  prescribed  tangential  stress  at  the  interface 
such  that  a  self  similar  problem  could  be  obtained.  He  then  employed  the  numerical 
solution  to  the  Boundary  Layer  Equation  in  an  inviscid  Orr-Sommerfeld  equation 
incorporating  the  effect  of  surface  tension,  and  a  stability  analysis  was  performed.  For 
the  case  of  constant  shear  stress,  his  analysis  revealed  the  existence  of  unstable  short 
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wavelength  surface  waves;  i.e.,  according  to  this  theory  the  abrupt  change  of  boundary 
conditions  at  the  nozzle  exit  would  indeed  be  responsible  for  the  atomization  of  the  jet. 

Hanratty,  T.  J.  (1983)  tried  to  solve  the  linear  momentum  equations  to  determine 
whether  small  amplitude  wavelike  disturbances  at  the  interface  will  grow  or  decay  when 
air  flows  concurrently  with  a  thin  liquid  film  in  an  enclosed  channel.  In  this  framework 
the  gas  flow  is  found  to  affect  the  stability  of  the  liquid  film  through  the  imposition  of 
pressure  and  shear  forces  at  the  interface.  He  described  that  the  prediction  of  the  pressure 
and  shear  stress  variation  along  a  wavy  surface  over  which  a  turbulent  gas  is  flowing  then 
becomes  the  central  problem  in  predicting  the  stability  of  the  liquid  film.  He  also 
described  that  a  slight  increase  of  the  gas  velocity  above  that  required  to  produce  two- 
dimensional  waves  causes  the  interface  to  break  into  a  pebbled  structure  with  a  distance 
between  crests. 

The  two-dimensional  waves  are  found. to  be  more  stable  on  liquids  of  larger 
viscosity  than  water;  i.e.,  they  exist  over  a  wider  range  of  gas  velocities.  He  described 
that  for  very  high  gas  velocities  wave-induced  pressure  variations  in  phase  with  the  wave 
height  can  be  so  large  that  they  overcome  the  stabilizing  influence  of  surface  tension. 
He  described  that  under  these  circumstances  small  wavelets  riding  on  top  of  the  roll 
waves  are  atomized  and  these  waves  are  generated  because  wave-induced  pressure 
variations  in  the  gas  that  are  in  phase  with  the  wave  slope  transfer  energy  from  the  air  to 
the  liquid  through  velocity  fluctuations  normal  to  the  interface.  He  described  that  an 
instability  is  predicted  when  the  rate  of  this  energy  input  exceeds  the  rate  of  dissipation  in 
the  liquid. 


However,  for  very  thin  film,  according  to  Hanratty,  wave-induced  velocity 
fluctuations  in  the  liquid  have  much  larger  components  in  the  tangential  direction  than  in 
the  normal  direction.  Consequently,  the  shear  stress  amplitude  can  be  an  important 
mechanism  for  transferring  energy  from  the  air  to  the  liquid.  The  capillary  ripples  appear 
on  thin  liquid  films  at  high  gas  velocities.  For  very  high  gas  velocities  for  all  film  flow 
rates,  the  destabilizing  influence  of  wave-induced  pressure  variations  is  found  to  be 
greater  than  the  stabilizing  influence  of  gravity  and  surface  tension;  very  rapid  wave 
growth  is  found  to  occur. 

Sattelmayer,  T.  and  Wittig,  S.  (1986)  showed  through  their  experiment  for  the 
airblast  film  disintegration  that  the  film  disintegration  becomes  more  uniform  with 
increasing  air  velocity  or  with  reduced  viscosity  and  surface  tension.  They  described  that 
the  disintegration  of  the  liquid  film  always  takes  place  close  to  the  atomizing  edge, 
regardless  of  the  liquid  properties.  They  described  that  small  liquid  ligaments  are  formed 
in  the  main  flow  direction  while  string-like  structures  are  observed  only  at  very  high 
viscosities  and  low  air  velocities.  At  high  liquid  viscosity,  the  air  is  found  to  be  unable  to 
atomize  the  liquid  immediately  after  separation  due  to  the  high  damping  capacity. 
Ligaments,  therefore,  are  formed  which  can  be  observed.  At  high  air  velocities,  the  liquid 
disintegration  is  found  to  take  place  in  a  very  small  zone  near  the  atomizing  edge. 
Immediately  after  separating  from  the  atomizing  edge,  the  liquid  is  forced  to  move  with  a 
strong  component  normal  to  the  main  flow  direction  of  the  air.  Subsequently,  the  liquid 
is  deflected  in  the  opposite  direction  by  the  aerodynamic  forces.  The  atomization  of  the 
liquid  begins  simultaneously. 
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Lin,  S.  P.  and  Kang,  D.  J.  (1987)  investigated  the  stability  of  an  axisymmetric 
liquid  jet  issued  into  an  ambient  gas  with  respect  to  spatially  growing  disturbances.  They 
assumed  that  the  liquid  is  an  incompressible  Newtonian  fluid  and  the  gas  is  assumed  to  be 
an  incompressible  inviscid  fluid.  After  the  corresponding  linearized  boundary  conditions 
are  applied,  they  obtained  normal  mode  solutions  which  contain  the  modified  Bessel 
functions.  From  their  mathematical  approach,  Lin  and  Kang  showed  that  the  onset  of 
atomization  is  due  to  the  pressure  fluctuation  which  causes  the  resonant  oscillation  of 
capillary  waves  at  the  liquid-gas  interface.  The  initial  disturbance  amplitude  at  the  nozzle 
exit  must  be  larger  than  a  threshold  amplitude  for  the  atomization  to  take  place 
immediately  at  the  nozzle  exit.  Otherwise,  atomization  may  also  be  caused  by  shear 
waves,  which  tend  to  produce  droplets  of  a  size  several  orders  larger  than  the  capillary 
length. 

Reitz  (1987)  first  extended  the  stability  analysis  of  the  atomization  of  fuel  jet 
injected  from  fuel  injector  in  internal  combustion  engine.  He  investigated  the  stability  of 
a  cylindrical  liquid  surface  to  perturbations  using  a  first  order  linear  theory.  He  imposed 
an  infinitesimal  axisymmetric  displacement  on  the  liquid  surface  and  ultimately  derived  a 
dispersion  equation  from  the  linearized  hydrodynamical  equations  for  the  liquid  and  gas. 

This  equation  relates  the  growth  rate  of  an  initial  perturbation  of  infinitesimal 
amplitude  to  its  wavelength.  The  relationship  also  includes  the  physical  and  dynamical 
parameters  of  the  liquid  jet  and  the  surrounding  gas  and  modified  Bessel  functions.  He 
predicted  that  the  presence  of  finely  atomized  drops  is  seen  at  the  spray  edge  near  the 
nozzle  in  high-  pressure  sprays  and  a  short  intact(unbroken)  region  of  the  jet  is  still  seen 
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close  to  the  nozzle  exit  in  high  pressure  sprays.  Later,  the  investigation  for  jet 
atomization  in  case  that  both  the  gas  phase  and  the  liquid  phase  velocities  are  non-zero 
had  been  done. 

Chuech,  S.  G.  and  Singhal,  A.  K.  (1991)  proposed  the  numerical  model  for 
primary  atomization.  They  formulated  for  the  one-dimensional  mean  flow  equations  for 
axisymmetric  jets.  They  assumed  that  the  gas  phase  flowfield  is  in  a  uniform  or  constant- 
velocity  profile  and  the  liquid  phase  flowfield  in  another  uniform  or  constant  velocity 
profile.  They  called  different  velocities  in  different  layers  as  step-varying  velocity 
profile.  They  derived  the  surface-wave  dispersion  equations  with  two  types  of  wave- 
models,  namely  linear  and  approximate  nonlinear  models.  From  the  calculation  with  the 
liquid  surface  and  jet  equations,  they  concluded  that  for  low  speed  jets,  only  the  large 
waves  grow  and  drops  are  formed  at  the  tail  of  the  jet;  and  for  high  speed  jets,  small 
waves  grow  at  the  liquid  gas  interface  and  provide  continuous  liquid  stripping.  They  also 
commented  that  the  wave-dispersion  equation  is  based  on  the  linear  small  perturbation 
theory  and  cannot  predict  jet  breakup  due  to  larger  (e.g.,  turbulence  generated) 
perturbations.  Further  study  of  liquid  breakup  criteria  is  needed. 

Lefebvre,  A.  H.  (1992)  experimented  with  jet  atomization  for  two-dimensional 
gas  phase  flow.  He  assumed  that  the  liquid  phase  flow  field  is  one-dimensional,  and,  the 
atomizing  air  impinges  on  the  liquid  sheet  at  small  angle.  This  method  of  air  impaction 
prevents  the  wavy-sheet  formation  that  in  classical  theories  of  liquid  sheet  disintegration 
is  generally  regarded  as  an  essential  first  step  in  the  atomization  process.  Before  any 
waves  have  time  to  develop,  the  liquid  sheet  emerging  from  the  nozzle  is  completely 
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disrupted  by  the  atomizing  air.  The  changes  in  the  angle  at  which  the  atomizing  air 
impacts  the  liquid  sheet  would  be  expected  to  affect  the  efficiency  of  the  atomization 
process. 

Samuelsen,  G.  S.  and  Stapper,  B.  E.  (1992)  experimented  with  the  breakup  of  a 
two-dimensional  liquid  sheet  in  the  presence  of  co-flow  air  with  emphasis  on  the  extent 
to  which  liquid  properties  affect  breakup.  The  liquid  (lower  surface  tension  compared  to 
water)  displays  smaller  scale  structures  on  the  sheet,  and  a  greater  tendency  to  breakup 
into  droplets  rather  than  form  streamwise  ligaments.  On  the  other  hand,  the  liquid  (higher 
viscosity  than  water)  stretches  farther  before  breakup,  which  results  in  the  streamwise 
ligaments  reaching  smaller  diameters.  The  breakup  time  of  the  sheet  is  dependent  on 
liquid  properties  and  relative  velocity  of  the  liquid  and  air.  The  liquid  (lower  surface 
tension  and  density)  demonstrates  decreased  breakup  time,  while  the  liquid  (higher 
viscosity  and  density)  exhibits  increased  breakup  time  relative  to  water. 

However,  the  physical  breakup  of  fuel  in  gas  turbine  combustors  actually  occurs 
three-dimensionally.  The  pressure  atomizer  and  airblast  atomizer  have  predominantly 
been  used  to  atomize  the  liquid  fuel  in  gas  turbine  engines. 

The  superior  combustion  and  reduced  emissions  of  the  airblast  fuel  injector 
relative  to  the  conventional  pressure  atomizer  is  well  established  by  Lefebvre  (1972). 
Unlike  pressure  atomizers  that  are  also  used  quite  extensively,  airblast  nozzles  feature  a 
significant  control  range  as  a  result  of  the  decoupling  of  atomization  quality  and  the  fuel 
mass  flow.  In  addition,  it  is  easier  to  control  the  spatial  distribution  of  droplets.  This 
superior  performance  has  led  to  the  widespread  adoption  of  the  airblast  fuel  injector  in 
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modern  high-pressure  ratio  gas  turbine  engines.  The  airblast  fuel  injector  produces  a  finer 
spray  because  of  the  thorough  mixing  of  the  fuel  and  air.  This  yields  a  flame  of  low 
luminosity  and  soot,  thus  resulting  in  relatively  cooler  liner  walls,  a  minimum  of  exhaust 
smoke,  and  lower  NOx  emissions.  Furthermore,  the  placement  of  the  fuel  droplets  in  the 
burning  zone  is  mainly  dictated  by  the  airflow  pattern  being  almost  insensitive  to  fuel 
flow  variations.  Thus  airblast  fuel  injectors  exhibit  lower  pattern  factors  and  therefore 
offer  an  important  practical  advantage  in  terms  of  turbine  nozzle  guide  vane  durability. 

Increasing  use  is  being  made  of  airblast  atomizers  in  stationary  combustion 
processes  and  for  fuel  preparation,  particularly  in  gas  turbine  combustion  chambers.  An 
almost  ideal  combination  results  in  cases  where  the  atomization  performance  can  be 
matched  to  the  required  operating  parameters.  These  are  mainly  the  mass  flow  provided 
for  the  combustion  at  a  given  load  point  and  the  pressure  loss  in  the  primary  air  needed  to 
stabilize  the  flame.  For  this,  it  is  advantageous  to  swirl  the  primary  air  used  for  flame 
stabilization  before  atomization  of  the  fuel. 

Aigner,  M.  and  Wittig,  S.  (1987)  experimented  with  the  atomization  of  the  liquid 
film  exposed  to  the  surrounding  air  coming  out  from  the  inner  airswirler  and  the  outer 
airswirler  in  gas  turbine  combustor.  They  observed  that  the  velocity  profiles  of  the 
atomization  air  and  the  properties  of  the  liquid  are  the  dominant  parameters  that 
determine  the  drop-size  distribution  generated  at  the  atomizer  edge.  They  also  observed 
that  the  shear  flow  at  the  nozzle  exit  depends  largely  on  the  swirl  and  counterswirl  of  the 
exiting  air  and  the  swirl  naturally  has  a  considerable  effect  on  the  flow  field  and  the 
droplet  path  downstream  from  the  nozzle.  As  explained  above,  the  sizes  of  the  droplets 
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produced  are  not  only  inversely  proportional  to  the  air  velocity  but  the  breakdown  is  also 
strongly  affected  by  the  flow  profile  of  the  atomization  air  at  the  atomization  edge.  An 
increase  in  the  swirl  velocity  produces  a  clear  reduction  in  axial  velocity. 

At  the  same  time,  the  pressure  loss  across  the  nozzle  increases.  With  an  identical 
pressure  loss,  however,  a  finer  atomization  would  have  been  attainable  only  by  utilizing 
counterswirl  because  the  high  additional  turbulence  intensity  is  generated  in  the  shear 
zone  downstream  of  the  atomization  edge.  The  volumetric  flow  and  the  viscosity  of  the 
liquid  do  not  show  a  direct  effect  on  the  droplet  size  produced  at  the  atomization  edge. 
With  a  proper  selection  of  swirl  elements,  and  particularly  with  counterswirl,  the  shear 
stresses  effective  on  the  atomization  edge  can  increase,  resulting  in  a  finer  atomization. 
Particularly  the  size  of  the  largest  droplet  classes  is  reduced.  In  all  cases,  the  amount  of 
air  required  by  the  atomizer  to  produce  a  constant  quality  of  atomization  is  significantly 
reduced  employing  swirl.  Swirling  the  atomization  air  in  order  to  stabilize  the  flame  can 
be  achieved  without  reducing  the  atomization  efficiency.  Flow  stall  at  the  atomization 
edge  should  be  avoided  by  reducing  the  swirl  angle. 

McDonell,  V.  G.  and  Samuelson,  S.  (1990)  utilized  a  two-component  phase 
Doppler  interferometry  to  investigate  the  two-phase  flow  produced  by  a  gas  turbine  air- 
blast  atomizer.  The  atomizer  is  operated  at  an  air-to-fuel  ratio  of  1 .0,  with  a  fuel  mass 
flow  of  0.00203  kg/s.  Methanol  was  employed  as  liquid  fuel  at  ambient  condition.  For 
the  relatively  low-pressure  drops  of  the  present  study,  the  air  flow  is  approximately 
incompressible.  The  fuel  is  injected  through  a  swirl  chamber  to  form  a  sheet  that  is 
immediately  exposed  to  air  on  both  sides.  The  outer  air  flows  through  a  swirler,  while  the 
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inner  air  passes  through  a  number  of  tangential  slots  to  the  central  passage.  The  helical 
vanes  of  55  degrees  for  swirling  are  used  to  atomize  the  liquid  fuel  and  form  the 
recirculation  zone.  The  complete  breakup  of  the  annular  fuel  sheet  into  drops  had 
occurred  at  the  atomization  edge. 

Rizk,  N.  K.  and  Mongia,  H.  C.  (1991)  used  data  from  McDonell,  V.  G.  and 
Samuelson,  S.  (1990)  to  verify  the  performance  of  a  numerical  model  for  the  initial 
atomization  in  a  gas  turbine  combustor.  That  model  was  applied  to  an  airblast  atomizer 
that  incorporated  a  short  prefilming  device  and  was  formulated  to  include  the  numerical 
simulation  of  the  airflow  through  and  exiting  from  the  atomizer  passages,  fuel  filming 
and  breakup.  According  to  Rizk  and  Mongia,  the  fundamental  processes  of  the  fuel 
atomization  include  liquid  formation  into  a  sheet,  its  subsequent  disintegration  into 
ligaments  and  drops.  The  prefilmer  geometry,  fuel  and  air  properties,  and  operating 
parameters  are  the  main  factors  controlling  the  film  thickness.  The  breakup  of  the  liquid 
sheet  into  ligaments  and  drops  involves  the  instability  of  the  liquid  stream  and  its 
tendency  to  form  waves  as  the  result  of  the  interaction  with  the  air.  Being  unstable,  the 
wave  amplitude  rapidly  increases  to  a  critical  value  at  which  the  liquid  disintegrates  into 
ligaments,  which  in  turn  break  up  into  drops. 

Rizk,  N.  K.  (1977)  shows  that,  with  an  increase  in  air  velocity,  the  liquid  sheet 
disintegrates  earlier  so  that  ligaments  are  formed  nearer  the  lip.  These  ligaments  tend  to 
be  thinner  and  shorter  and  disintegrate  into  smaller  drops.  The  balance  of  forces  caused 
by  gas  pressure,  surface  tension,  liquid  inertia,  and  viscosity  controls  the  breakup 
mechanism  of  the  sheet.  It  is  shown  that  the  breakup  of  planar  liquid  sheets  due  to  the 
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wave  growth  concept  occurs  as  the  total  growth  of  the  wave,  which  is  given  by  the  natural 
logarithm  of  the  ratio  of  wave  amplitude  to  amplitude  of  initial  disturbance. 

However,  there  is  no  evidence  of  ligaments  existing  in  the  photograph  of  the  spray 
structure  by  McDonell  and  Samuelsen  (1990). 

The  mathematical  investigation  for  annular  layers  with  density  variation  and 
swirling  flows  has  been  made  by  a  few  researchers.  Reshotko,  E  and  Monnin,  C.  F. 
(1965)  studied  the  stability  of  an  incompressible  two-fluid  wheel  flow  to  infinitesimal 
helical  disturbances.  In  their  analysis,  Reshotko  and  Monnin  found  that  as  the  density 
ratio  of  inner  to  outer  fluid  increases,  the  calculated  growth  rates  increase  and  approach  a 
limiting  value  as  the  density  ratio  becomes  infinite  and  growth  rates  also  increase  with 
increasing  axial  as  well  as  azimuthal  wave  number. 

Fung,  Y.  T.  (1974)  investigated  the  instability  growth  rate  for  three  regions  of 
fluids  with  the  cylindrical  interfaces  in  a  constant  angular  and  axial  velocity  field  and  step 
type  density  distribution.  He  derived  the  dispersion  equations  by  asymptotic  form  of 
Bessel  functions  by  matching  the  kinematic  and  dynamic  boundary  conditions  at  both  the 
interfaces.  He  observed  that  the  growth  rate  increases  when  the  density  in  middle  region 
is  larger  than  those  in  other  regions  and  those  in  other  regions  are  the  same. 

Fung,  Y.  T.  and  Kurzweg,  U.  H.  (1975)  investigated  the  inviscid  instability  of 
heterogeneous  swirling  flows  with  radial-dependent  density  and  obtained  the  secular 
relations  for  the  instability  growth  rates  for  several  different  flow  configurations.  They 
observed  that  very  stable  flows  occur  when  the  angular  and  axial  velocity  components 
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have  no  radial  variation  and  simultaneously  the  density  increases  with  radius  as  is  the 
case  in  a  typical  centrifuge. 

Fung,  Y.  T.  (1983)  also  derived  the  mathematical  formulation  for  rotating  flow  for 
three  annular  layers  with  density  variation.  Only  azimuthal  modes  are  considered. 
Angular  velocity  in  a  layer  is  not  same  as  that  in  another  layer. 

Meyer,  J.  and  Weihs,  D.  (1987)  investigated  the  capillary  instability  of  an  annular 
liquid  sheet.  In  their  analysis,  the  researchers  suggested  that  the  annular  jet  behaves  like  a 
two-dimensional  liquid  sheet.  The  most  unstable  perturbations  are  antisymmetric  and 
their  growth  rate  increases  as  the  jet  thickness  decreases.  Therefore,  an  annular  liquid  jet 
with  a  sufficiently  small  ring  thickness  will  disintegrate  into  spherical  shells  much  faster 
than  a  full  liquid  jet  disintegrates  into  drops. 

Lin,  S.  P.  and  Zhou,  Z.  W.  (1992)  investigated  the  atomization  of  a  liquid  jet 
emanating  from  a  hole  nozzle  into  the  ambient  gas  for  the  application  to  internal 
combustion  engines.  They  used  the  linear  stability  analysis.  They  assumed  that  the  fluid 
is  inviscid,  incompressible  with  surface  tension.  Only  axial  velocity  components  were 
considered  in  their  analysis.  They  showed  the  spatial  growth  of  the  disturbance  amplitude 
in  the  flow  direction  when  the  instability  occurred.  They  proposed  that  the  relation 
between  the  radius  of  droplet  a  and  the  wavelength  of  disturbance  X  is  nR2X  =  4na3 1 3 
where  R  is  the  radius  of  a  circular  cylindrical  nozzle.  They  showed  that  the  surface 
tension  has  an  impact  on  damping  of  growth  of  the  disturbance  amplitude. 

In  summary,  the  fuel  sheet  moves  axially  and  circumferentially  against  the 
surrounding  gas  which  is  moving  axially  and  tangentially  with  a  different  speed  in  the  gas 
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turbine  combustor.  That  means  that  the  fuel  sheet  moves  helically  just  after  the  tip  nozzle 
in  the  gas  turbine  combustor.  To  investigate  the  instability  of  the  fuel  sheet,  both  the 
instability  of  the  axial  flow  annular  jet  and  the  instability  of  the  rotating  flow  annular  jet 
must  be  investigated  in  advance  in  order  to  investigate  and  understand  the  helical  flow 
annular  jet.  To  understand  the  physical  phenomena  of  instabilities,  the  investigation  of 
Meyer,  J.  and  Weihs,  D.  (1987)  is  helpful  for  the  axial  flow  annular  jet,  and  the  Reshotko 
and  Monnin  (1965)  investigation  is  helpful  for  the  rotating  flow  annular  jet.  The 
investigation  of  Kurzweg  (1970)  is  especially  helpful  in  understanding  the  physical 
phenomena  for  the  helical  mode  of  instability  in  swirling  flow  between  concentric 
cylinders.  Also,  the  investigation  of  Fung  (1983)  is  helpful  in  understanding  the  physical 
phenomena  for  the  instability  of  three  rotating  annular  layers.  The  extension  of  linear 
instability  analysis  by  Reitz  (1987)  and  Lin,  S.  P.  and  Zhou,  Z.  W.  (1992)  into  the 
application  of  internal  combustion  engine  has  brought  the  idea  of  the  extension  of  linear 
instability  analysis  in  gas  turbine  combustors. 


CHAPTER  2 
MATHEMATICAL  FORMULATION 


The  front  view  of  the  fuel  atomizer  is  shown  in  figure  2  and  side  view  of  the  fuel 
atomizer  is  shown  in  figure  3.  The  fuel  leaves  the  nozzle  as  an  approximate  annular  sheet 
surrounded  inside  and  out  by  air.  The  schematic  of  the  spray  structure  in  the  experiment 
is  shown  in  figure  4. 

The  fuel  injection  into  gas  turbine  combustor  was  modelled  in  this  dissertation  as 
a  helical  disturbance  of  the  annular  liquid  fuel  sheet  exposed  to  the  air.  The  fluids  are 
assumed  incompressible  and  inviscid.  Wheel  flow  is  applied  to  the  air  and  the  liquid  fuel. 
However,  each  has  a  different  angular  velocity.  The  axial  velocities  of  the  inner  air, 
liquid  fuel  and  outer  air  are  different  from  one  another. 

The  experimental  data  used  in  mathematical  formulation  are  provided  by 
McDonell  and  Samuelsen  (1990)  at  the  University  of  California,  Irvine  Combustion 
Laboratory.  From  the  experimental  data,  the  following  data  for  investigating  the 
instability  of  an  annular  fuel  sheet  are  calculated: 

1 .  axial  velocity  components  of  both  the  surrounding  gas  and  fuel  sheet  at  atomizer  exit 

2.  swirling  velocity  components  of  both  the  surrounding  gas  and  fuel  sheet  at  atomizer 
exit 

3.  fuel  sheet  thickness 
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4.  density  of  fuel,  density  of  air 

5.  surface  tension 

6.  inner  radius  and  outer  radius  at  interfaces  between  the  annular  fuel  sheet  and 
surrounding  gas 

To  understand  the  behavior  of  the  instability  mechanism  in  swirling  flows,  it  is 
necessary  to  consider  the  mathematics  of  the  Navier-Stokes  equations,  the  continuity 
equation,  and  incompressibility  and  solve  these  governing  differential  equations.  The 
onset  of  instability  can  be  determined  by  assuming  temporal  or  spatial  disturbances  whose 
behavior  is  governed  by  the  linearized  forms  of  the  governing  nonlinear  equations.  It  has 
been  accepted  that  the  method  of  normal  modes  is  complete  and  well  explains  the  linear 
wave  behavior  in  fluid  mechanics.  The  linearized  partial  differential  equations  under 
consideration  are  cyclic  with  respect  to  time  and  to  two  spatial  variables.  This  strongly 
suggests  they  be  solved  by  Laplace  (time)  and  Fourier  (space)  transform  methods.  The 
differential  equations  after  transformation  are  ordinary  and  in  general  complex  both  with 
respect  to  the  temporal  and  spatial  parameters. 

The  solution  to  the  governing  equations  has  as  a  complex  eigen  growth  rate  co  = 
cox  +  ico  j.  Whenever  co,  is  positive,  instability  occurs.  The  breakup  of  the  annular  fuel 
sheet  is  caused  by  the  increase  of  growth  rate  ( co ,).  We  wish  to  derive  the  instability 
mechanism  at  the  interface  between  the  surrounding  air  and  the  annular  fuel  sheet 
exposed  to  the  surrounding  air.  Both  the  surrounding  air  and  the  annular  fuel  sheet  are 
swirling  and  proceeding  downstream  just  after  the  nozzle.  Therefore,  they  have  the 
steady-state  radially  dependent  velocity  field  V[0,rQ(r),W(r)].     The  fluid  under 
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consideration  is  assumed  inviscid,  incompressible  and  non-heat-conducting  but  has  a 
radially  dependent  density  po(r).  Euler's  equations  of  motion,  the  incompressibility 
condition  and  the  continuity  equation  governing  such  a  swirling  flow  in  cylindrical 
coordinates  (r,  9,  z)  are,  respectively, 
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To  derive  the  instability  of  two  phase  flow,  we  postulate  an  initial  set  of  solutions 
based  on  the  summation  of  the  steady-state  flow  and  infinitesimal  time-dependent 
fluctuations,  namely, 

Vr  =  u(r,0,z;t) 

Ve  =  rQ(r)  +  v(r,0,z;t) 


Vz  =  W(r)  +  w(r,e,z;t) 


(2.4) 


P  =  po(r)  +  p(r,0,z;t) 


P  =  Po(r)  +  p(r,0,z;t) 

Substituting  equations  (2.4)  into  equations  (2.1),  (2.2)  and  (2.3),  neglecting  the 
gravitational  effects,  and  following  the  normal  linearization  procedure  of  dropping  second 
and  higher  order  terms  in  the  assumed  infinitesimal  perturbations  lead  to  the  linearized 
differential  equations 


du    u    1  dv    dw  . 

—  +  -  +  +  —  =  0 

dr    r    r  30  dz 


(2.5a) 


d.p  +  (DPo)u  =  0 


(2.5b) 


dr 


+  pQ{d.u-2Q.v)-prQ.2  =0 


(2.5c) 


^  +  rpo[d.v  +  D\rQ)u]  =  0 


69 


(2.5d) 
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p0[d.w  +  (DW)u]  =  0 


(2.5e) 


where 


D  =  —,D'  =  D  +  -    and    d.=  —  + 


dr  r  dt 


d_ 
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To  reduce  the  number  of  perturbation  variables,  it  is  desirable  to  eliminate  v,  w 
and  p  because  u  and  p  are  primarily  related  to  boundary  conditions  and  interfacial 
conditions. 

To  eliminate  v,w  and  p  from  eq.  (2.5a)  to  (2.5e),  we  convert  eq.  (2.5c)  into  the 
terms  including  d.  ,  then  we  obtain 


d,[p0(d.u  -  2Qv)]  -  d.(prQ2)  =  0 


or 


d.~  +  p0d.2u - 2pQd.v  - rQrd.p  =  0 . 


dt 


(2.5c.l) 


We  multiply  —  into  eq.  (2.5d),  then  we  obtain 
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^  fU  2Qp0[d.v]  +  p0  [4Q2  +  2Qr  ^]fi  =  0 . 
r  dG  dr 


(2.5(1.1) 


From  eq.  (2.5b),  d.p  =  -(Dp0)u  is  derived  and  this  equation  is  substituted  into 
(2.5c.  1)  and  we  obtain 


rln 

d.  —  +  p0d,2u  -  2pQd,v  +  rQ2(Dp0)u  =  0 .  Once  this  equation  is  added  into  (2.5d.l), 
dr 


we  obtain 


,  dp  2Q  dp  ,2.  A 
a,  —  +  —  +  pnd,  u  +  Ow  =  0 


(2.6) 


where     0  =  —  [rQ2 +  4/?0Q2  +  2/y*2—  ] 
p0         i/r  dr 


We  differentiate  eq.(2.5d)  with  respect  to  #  and  multiply  \,  then  obtain 
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We  differentiate  eq.(2.5e)  with  respect  to  z  ,  then  obtain 


?±+Ptld.*(DW&=0 

oz  oz  oz 


(2.5e.l) 


Once  (2.5d.2)  and  (2.5e.l)  are  added,  then  we  obtain 
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1  & 'p  d2  p  r  ,  A  oO  dw.  D' .  ^du  /T.TI,^du1  . 
r  d&      az  r  06    dz       r         dO  dz 


(2.5f.l) 


From  eq.  (2.5a),  we  obtain 


1  dv  ^dw  _  du  u 
r  d9    dz       dr  r 


\_d_ 
r  dr 


(ru) 


Adding  eq.  (2.5f.l)  to  eq.  (2.5f.2),  we  obtain 


1  d2p    d2 p       .   ,  .  1  d  .     ,    D'  .  ^  du    .       du^  n 
+  ^  +  P^rd.{-—(ru)}  +  — (rQ)—  +  (DW)—\  =  0 
r  dd      dz~  r  or  r         o9  dz 


(2.5f.2) 


(2.7) 


The  coefficients  from  eq.  (2.5a)  to  eq.(2.7)  depend  only  on  the  radial  coordinate 
r,  while  the  remaining  independent  variables  0,  z  and  t  enter  only  in  the  partial 
derivatives. 

This  allows  one  to  assume  that  (2.5),-  (2.6)  and  (2.7)  admit  solutions  having 
an  exponential  behavior  in  z,9  and  t.  This  assumption  reduces  all  these  partial 
differential  equations  to  ordinary  differential  equations.  Following  the  normal  analysis 
indicated  in  Chapter  1,  we  apply  the  Laplace  transform  with  respect  to  time  t  and  the 
Fourier  transform  with  respect  to  the  spatial  variables  (z,0).  The  transforms  have  the 
form 


{u,p}(k,m,co;r)  =  j  \  ]{u,p}(r,0,z;t)el{la+m9-a,)ciedzdt 


(2.8) 


0  -« 0 
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which  is  equivalent  to  the  standard  transformation  of  assuming 


{p,u}  =  {p(r),u(r)}e 


i(kz+m9-w() 


We  are  here  considering  only  temporal  disturbances  so  that  the  axial  wave  number  k  is 
real  while  the  amplification  factors  =  cor  +  icOj  is  complex;  m  represents  the  azimuthal 
wave  number  and  takes  on  only  integer  values  because  of  the  periodicity  of  the  solution 
in  the  9  direction.  Using  transform  (2.8)  in  equations  (2.6)  and  (2.7)  leads  to  the  two 
ordinary  differential  equations 


NDp  + 


2mQ. 
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where 
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39       dz        dr  r 


and 


TV  =  kW  +  mQ.  -  co 


CHAPTER  3 
BOUNDARY  CONDITIONS 

3.1    Inlet  Conditions 

The  boundary  conditions  for  the  airblast  atomizer  delivering  swirling  airflows  on 
either  side  of  an  annular  fuel  sheet  are  described  in  this  chapter.  Inner  and  outer  counter 
clockwise  swirled  airstreams  atomize  and  distribute  fuel  of  a  thin  swirling  annular  sheet. 
Since  the  effect  of  the  viscosity  is  regarded  negligible,  the  flow  field  consists  of  three 
layers  as  seen  in  figure  2.  The  layers  have  discontinuous  densities  and  velocities  since  the 
effect  of  the  viscosity  is  regarded  as  negligible.  The  inner  layer  air  stream  core  includes 
the  axis.  Next  is  the  fuel  sheet  layer.  Beyond  the  fuel  sheet  is  the  outer  air  stream  layer. 

The  nozzle  tip  imparted  a  radial-inflow  component  to  the  outer  swirled  stream  by 
forming  an  angle  with  the  axis  of  the  nozzle.  The  front  view  of  the  nozzle  section 
surrounded  by  the  air  is  shown  in  figure  2.  The  side  view  of  the  atomizer  head  section  is 
shown  in  figure  3.  The  effect  of  the  certain  angle  with  which  the  outer  swirled  stream 
intersects  the  annular  fuel  sheet  is  not  considered  in  this  dissertation.  Therefore,  it  is 
assumed  that  the  nozzle  tip  is  parallel  to  the  axis  of  the  nozzle. 

The  mean  axial  velocities  and  circumferential  velocities  are  calculated  in 
accordance  with  the  dimensions  and  configurations  of  the  airblast  atomizer.  The  data 
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from  the  experimental  setup  is  used  to  calculate  the  axial  velocities  of  both  airstream  and 
fuel  sheet.  The  pressure  drop  for  both  fuel  AP^c/  and  inner  air  A/V.  given  from  the 

experimental  setup  are  used  to  calculate  the  axial  velocities  of  both  airstream  and  fuel 
sheet.      From  the   Bernoulli  equation,    A/^w  =  —  p^W^,2 ,   the   fuel  velocity, 


W  = 

"fuel 


\2AR~,  , 

 —  is  obtained.  pM  is  known.  The  helical  vane  angle  of  55    is  used  to 


y 

determine  the  fuel  swirl  VtlJM)  =    Wfuel*  tan(55°).   Qfuel  =      tl)  is  obtained  where 


'fuel 


rfuel  is  the  radial  distance  from  the  axis  to  the  centerline  of  the  fuel  layer. 


Similarly,  Vt 


2AP.  . 


air(inner) 


is  obtained  and  used  to  calculate  VMn . .  and 


{inner) 


Vair(inner)  =  9(inner)  +        dinner)    and  V ^   =      Wgir  *  tan(55°)  due  tO  helical 


vane  angle  of  55  Wair  =  V^r 


Vl  +  tan(550)2  ' 


^air  ~  — =   g(""'er>  ,  where  rmner  and  rouler  are  the  radial  distances  from  the 


r  r 

inner  outer 


axis  to  the  centerline  of  the  inner  air  layer  and  the  outer  air  layer,  respectively.  Vg( 


outer) 


is  obtained  from  V. 


d(outer) 


—  Vevnnery   Kr{outer)  *  calculated  from  IseSSL .  ^is 


tan(55u) 


approximately  regarded  as  10  times  density  at  atmospheric  pressure  in  gas  turbine 
combustor  inlet.  pfuel  is  the  standard  value  at  room  temperature. 
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In  summary,  axial  velocity  components  {Wair{imer),W ^,1?^^)  and  angular 
velocity  components  {0.air,D.fuet)  and  (pair,  pfuel)  are  now  setup.  The  values  are  as 
follows: 

A,,«=0-0122  Slcm' 

^r)=1.5010E+03c*/5 

n^^-22.8664E+03  Is 

0.787  g/cm3 
WMI=  1.4567E+03  cm  Is 
QMI  =  9.5649E+03  Is 

PaHoulerr  °-0122  Slcm' 
^,w=5.0836E+03  cmls 
Qa,r(OUIer)=22.8664E+03  Is 

In  this  study,   Wfiiel  >  Wair{ouler)  and  Wmr(mner)  is  only  considered  since  the  known 

solutions  which  have  W2  >  W3  and  Wi  are  available  to  investigate  the  validation  of  the 

computational  results  of  the  mathematical  development,   fi^  ,ow.  is  set  equal  to  zero 

since  the  radius  is  much  longer  to  the  combustor  wall  than  to  the  outer  interface  between 
the  annular  fuel  sheet  and  the  surrounding  outer  air  and  the  region  affected  by  ft^ .  is 


is  negligible  in  comparison  with  the  region  which  is  not  affected  by  0.air(outer)  in 
combustor. 

3.2    Interfacial  Conditions 

u ,  the  radial  perturbation  velocity,  is  finite  at  both  the  axis  and  the  outerwall  of 
the  combustor. 

The  use  of  non-smoothly  varying  profiles  can  be  sometimes  justified  as  an 
approximation  to  related  smoothly  varying  ones.  The  analysis  of  the  instability  for 
discontinuous  velocity  and  density  distributions  allows  simple  solutions  to  the  governing 
differential  equation  in  the  individual  regions  which  each  possess  continuous  simple 
varying  velocity  and  density  distributions. 

3.2.1    Kinetic  Interfacial  Condition 

The  fluid  regions  are  three  circular  rings  or  regions;  inner  air  region,  annular  fuel 
sheet  region  and  outer  air  region.  Consider  region  1  (inner  air)  and  region  2  (annular 
fuel  sheet)  with  an  interface  at  r  =  R.  Integrating  (2.10)  from  R  -  £  to  R  +e  in  the 
adjoining  subintervals  of  the  interface  at  r  =  R  and  assuming  that  the  mean  value  exists 
for  all  the  quantities  of  the  flow  and  that  their  derivatives  should  be  bounded  at  the 
interface.  We  obtain 
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R+e 


\ip0[k(DW) 


mD'(rQ) 


-ND']udr  =  0 


R-e 


(3.2.1) 


where 


e  - 


thickness  of  sheet 
R 


When  eq.  (3.2.1)  is  divided  by  N  ,  we  obtain 


R+£ 


\i^[k(DW)  +  ^-^-ND']udr  =  0. 


R-e 


N 


(3.2.2) 


When  eq.  (3.2.2)  is  divided  by  ip0,  we  obtain 


R+e  ,-.*/•  R-e 

—  [£(ZW)  +  * — ^]d>-    £»w  =  0. 


Here, 


r  w  r2wQ1  ,  ,    f  w  , 

—  [  Jar  =  0,  and     —  dr  =  0  are  consistent 

BJ  N     r  J  r 

R-e  R+e 


Therefore, 


R+e  R+e 

lim^  \R  \[m(DQ)  +  k(Dw)]dr-  {(Du)dr  =  0  or 

~°    N       R-e  rU 


N 


[m(Q2  -  Q, )  +  k(W2  -  Wx )]  -  [u2  -  ux  ]  =  0 


(3.2.3) 


where  — 
N 


*    2  A,    N2     N  *  N,-N 


u2  -  «, 


is  also  derived  from  eq.  (3.1). 
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Ultimately,  -(—  +  —) 


2  N,  AT, 


M2  ""I 


leads  to 


as 


-  =0 


is  defined 


(3.2.4) 


which  is  a  jump  condition  at  the  interface.  The  quantity  j^j  >  physically  interpreted  as  a 

perturbation  displacement,  must  be  nowhere  discontinuous  at  the  interface  or  anywhere 
inside  domain  of  the  flow.  Condition  (3.2.4)  is  also  obtainable  by  the  following  physical 
argument: 

Let  the  radial  perturbation  displacement  be  r\ .  The  analysis  of  the  normal  mode 
method  allows  the  perturbation  quantities  to  admit  solutions  with  periodic  behavior  in  z, 
6  and  t  and  with  radially  dependent  amplitudes.  Hence  we  assume 


\(kz+m8-w() 


fj(r,0,z;t)  =  T](r)e 


According  to  the  chain  rule 


(3.2.5) 


drj  _  dr]  ci]  dr  cr]  dQ  ci]  dz 
dt~  a     or  dt    39  dt     &  dt 


dt] 
~~dt 


=  i(-co  +  mQ.  +  kW)r)(r)ei{la+m0-a*) 


^-  =  iNn(r)eiUa+ma-M) 
dt  A 


(3.2.6) 


By  definition,  the  total  derivative  of  the  radial  perturbation  displacement  is  equal  to  the 
radial  perturbation  velocity  which  has  the  form 
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u  =  u(r)e 


(3.2.7) 


Equating  (3.2.6)  and  (3.2.7),  one  finds, 


u       .  u 


(3.2.8) 


Equation  (3.2.4)  then  follows  from  the  physical  argument  that  the  normal  component  of 
the  perturbation  displacement  must  be  continuous  across  the  material  fluid  interface 
inside  the  flow  domain. 

3.2.2  Dynamic  Interfacial  Condition 

The  general  dynamic  interfacial  condition  can  be  derived,  using  the  integration 
technique  across  the  interfaces  in  the  similar  way  used  for  kinematic  interfacial 
condition.  When  equation  (2.9)  is  divided  by  N ,  it  can  be  rearranged  to  the  form 


(3.2.9) 


Applying  the  integration  technique  across  the  interface  at  r  =  R,  we  obtain 


R-e 


(3.2.10) 


Rearranging  eq.  (3.2.10),  that  becomes 


When  integrated  becomes 


(3.2.11) 
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3.2.3  Surface  Tension 


When  we  consider  surface  tension  a j  acting  at  the  interface  j  located  between 

two  rings  of  fluids  at  r  =  Rj.  The  equation  of  motion  for  the  mean  flow  in  the  radial 
direction  is 


dR 


dr 


'  1 

--/V"2=-l7}(-)^-*,) 


Rj 


(3.2.12) 


where  S(r  -  AV )  is  a  Delta  function,  R    is  the  position  as  well  as  the  principal  radius  of 

curvature  for  the  disturbed  interface. 

Once  the  flow  is  disturbed,  the  principal  radii  of  curvature  on  r-  6  and  r-z  planes 
are  specified  as 


1       1      1     1  c?r  &r 

+  ■ 


Rrg    Rn    r    r2  dd1  dz1 


(3.2.13) 


and  the  linearized  perturbation  equation  of  motion  in  the  radial  direction  yields 


^  +  f  +  Po(dM-2Qv)-(p0+p)rQ2=-fjTJ(^-  +  ^  )5{r-Rj)  (3.2.14) 


Subtracting  (3.2.12)  from  (3.2.14),  one  finds 


f  +  p0(d.u-  2Qv)  -  prQ2  =  £  Tj  [-L  -  ( J-  + 1  )]S(r  -  R; ) 

W  J=I  A,  K,a  Krz 


(3.2.15) 
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In  the  absence  of  surface  tension,  equation  (3.2.15)  becomes  identically  equal  to  equation 
(2.5c).  Following  the  normal  mode  analysis  and  assuming  infinitesimal  disturbances, 
we  let 


r  =  Rj  +  ?,=  RJ  +  TKRJ)ei 


(kz+mO-cut) 


(3.2.16) 


at  the  interface  j.  Equation  (3.2.13)  can  be  linearized  by  the  rearrangement  of  each  term 
in  that  equation,  we  simplify 


1  1 


r     R  +  fj  R 


3r 


since  R»  tj.  With  the  following  expressions, 


—  =  imrj{Ri)e,(la+m9-a') 
d6  1 

~  =  -k2T1{Ri)e,^m9-^ 


The  bracket  in  eq.  (3.2.15)  becomes 


1,1  1x1  1 


i(kz+m6-aA) 


(3.2.17) 


(3.2.18) 


Here    rjR  )e'  '+m       can  be  expressed  as  -i  —    .  When  eq.  (3.2.15)  is  transformed 

\NJj 

into  ordinary  differential  equation  and  multiplied  by  N ,     eq.  (3.2.15)  becomes 

NDp+        p  =  m  _  PgN- )u_NYj-L.(l-m>-  Kj2  i^rlKr-R)}         (3 .2. 1 9) 

/-i  Rj  \Ns j 
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Dividing  eq.  (3.2.19)  by  N  and  applying  the  integration  technique  across  the  interface  at 
r  =  R ,  the  first  term  in  right  hand  side  leads  to 


R+e 


\i{D(p0m2)}  [±)dr 


R+c  ,  n 

since     ji{p0(3Q2 -N2)}  \-^-\dr  =  0.  Ultimately  eq.  (3.2.19)  becomes 

R-s 

[[Dp - i{D(pQrQ2 )-^(\-m2-K2)[^\  S(r-  R)}]dr  =  0 


or 


(3.2.20) 


CHAPTER  4 
ANALYTICAL  SOLUTIONS 

In  Chapter  4,  dispersion  relationships  are  derived  for  both  Reshotko's  two  wheel 
flows  and  for  annular  swirling  flows  with  discontinuous  velocity  and  density  profiles 
subject  to  arbitrary  infinitesimal  perturbations.  To  derive  the  dispersion  relationship,  it  is 
first  necessary  to  solve  the  system  of  equations  consisting  of  eq.  (2.9)  and  eq.  (2.10)  with 

the  boundary  conditions.  The  type  of  dispersion  relationship  is  F(^,m,kR,^)  =  0. 

The  growth  rate  co  depends  on  m  and  k .  The  growth  rate  (<».)  indicates  the  strength  of 
the  instability,  m  is  the  circumferential  wave  number  around  the  circular  interface 
disturbed,  while  k  is  the  axial  wave  number  around  the  axisymmetric  interface  disturbed. 
The  schematic  diagram  showing  surface  waves  on  an  annular  sheet  and  a  full  liquid  jet  is 
shown  in  figure  5.  When  a  helical  disturbance  occurs  on  the  infinite  cylindrical  core  jet, 
the  axisymmetric  disturbance  wave  occurs  in  the  axial  direction  and  the  tangential  waves 
anti-symmetrically  occur  in  the  tangential  direction  at  the  same  time.  Similar  phenomena 
occurs  on  the  annular  fuel  sheet  when  a  helical  disturbance  occurs  on  the  infinite  annular 
sheet. 

Consider  one  ring  of  fluid  with  constant  density  pJt  constant  angular  velocity  Q. 
and  constant  axial  velocity  Wj  in  region  j.  O  becomes  4Q2.  Eq.  (2.9)  and  (2.10)  are 
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side  view  (k  mode) 


front  view  (m  mode) 


fully  core  jet 


side  view  (k  mode) 


front  view  (m  mode) 
annular  sheet 


.  Schematic  Diagram  showing  surface  waves  on  an  annular  sheet 
and  a  fully  core  jet 


repeatedly  written  as 


2mQ  -, 

NDp  +  p  =  i(p0Q>  -  p0N-)u 

r 


(4.1) 


i—  +  k-)p  =  ip0[{k(DW)  +  i — J-}u-N(D  U)] 

r  r 


(4.2) 


To  eliminate  all  other  quantities  in  favor  of  the  pressure  disturbance  amplitude, 
eq.  (4.1)  can  be  reduced  to 


iPjQj\4  -  nf  )u  =  QjirijDp  +  -yp)  or 


i  ItYl 

ipjQj(4  -  nf)u  =  (rijDp  +  —  p) 


Eq.  (4.2)  can  then  be  reduced  to 


(4.3) 


(^T  +  k2)P  =  ip0[^u-N(D'u)]  or  &  +  k2)p  =  iponj[^-nJD-)u 
r  r  r  r 


From  (4.3),  u  is  rearranged  as 


-in. 


u  = 


Dp- 


2mi 


(4.4) 


(4.5) 


L)  u  =  —  +  -  m  eq.  (4.4)  is  expressed  as 
dr  r 


Du  = 


— — j-  -(—  Dp  +  njD'p  +  ^-Dp) 

pp.j{A-nj  )    r  r 


(4.6) 


Therefore,  eq.  (4.4)  is  rearranged  as 
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{k2  +  r^-)p  =  — —^—r  D  p  — -J-^r~  d'p  +  p  or 


D>p  +  ^{k\^)  +  ^}P  =  0 
r  n :  r- 


n,  -  4  , 
Here  —  is  defined  as  q ,  .  Once  eq.  (4.7)  is  multiplied  by  r2 ,  it  becomes 


(4.7) 


r2D2p  +  rDp-(k2q/r2  +m2)p  =  0 

It  is  called  the  disturbance  equation. 
Eq.  (4.8)  corresponds  to  modified  Bessel  function,  where 


2-  2   2   ,  2- 


P  =  ciIm(kqjr)  +  c2Km(kqjr) 


(4.8) 


where  q}  - 


(4.9) 


Alternatively,  eq.  (4.8)  corresponds  to  any  linear  combination  of  Bessel  functions  of 
the  first  and  second  kind  of  the  proper  argument,  ikq,r  where 


P  =  eJm(ikqJr)  +  jHmm(kqjr) 


where  qj  = 


(4.10) 


n. 


Although  eq.  (4.10)  is  used  less  as  a  general  solution  than  eq.  (4.9),  it  was  used  by 
Reshotko  (1965)  for  deriving  the  mathematical  theory  of  two-wheel  flows  with  different 
densities.  Reshotko  (1965)'s  mathematical  theory  will  be  used  in  Chapter  5  to  validate 
the  numerical  code.  Therefore,  the  mathematical  theory  with  eq.  (4.10)  for  two-wheel 
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flows  is  derived.  However,  eq.  (4.9)  is  preferably  used  to  derive  the  mathematical 
expressions  for  the  instability  of  the  annular  sheet  exposed  to  the  surrounding  gas. 

4.1  Reshotko's  Solutions  for  Two- Wheel  Flows 

As  shown  in  Figure  6,  the  wheel  flow  consisted  of  two  immiscible  and  rotating 
fluids  separated  by  a  cylindrical  interface.  The  flow  was  considered  to  be  only  azimuthal 
and  solely  a  function  of  the  radial  coordinate.  There  is  no  axial  velocity  component  and 
annular  velocity  in  the  inner  fluid,  Q,  is  equal  to  the  angular  velocity  in  the  outer  fluid, 
fi, .  And,  px  in  the  inner  core  liquid  is  larger  than  p2  in  the  outer  liquid.  In  the  inner 
liquid,  the  second  term  from  eq.  (4.10)  is  inadmissible  so  that 

Px=exJm{iqkr)  (4.11) 
In  the  outer  fluid,  the  first  term  from  eq.  (4.10)  is  inadmissible  so  that 
P2=f2Hm(l)(iqkr)  (4.12) 
e,  and  f2  are  unknowns.  Two  conditions  are  required  to  solve  for  those  unknowns.  The 
remaining  conditions  are  that  the  radial  displacement  must  be  continuous  at  the  interface 
and  that  the  normal  stress  must  be  continuous  at  the  interface.   These  conditions  are 
treated  consecutively.  Eq.  (3.2.4)  for  kinetic  interfacial  condition  and  eq.  (3.2.11)  for 
dynamic  interfacial  condition  are  applied  to  the  interface  between  the  inner  core  liquid 
and  outer  fluid.  To  the  general  solution  for  u ,  eq.  (4.5)  is  used.  Eq.  (4.5)  includes  Dp 
and  p.  When  the  derivative  term  of  eq.  (4.1 0)  with  respect  to  r  and  eq.  (4.10)  are 
substituted  into  eq.  (4.5),  eq.  (4.5)  becomes 
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Figure  6.  Geometry  of  two-fluid  wheel  flow 
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t,0*      sr'/-  ,  x    2mA J, (iqkr)^ 


-r.cy      .  TT  '.  ,  .  2mR,H(iqkr)^ 

f  [(- " m)Hm  (iqkr)  1   mK  *  J) 

£2  r 


(4.13) 


The  interfacial  conditions  for  eq.  (3.2.4)  and  eq.  (3.2.1 1)  become 


CO 


*i  [(-  "  m)Jm  {iqkRx  )-2mJm  {iqkR, )]  - 


CO 


f2  [(-  -  m)Hm  (iqkR,  )-2mHm  (iqkRx )]  =  0 


and 


i  CO  ' 


{co-m)  Q 


Pi. 


[  Mi~m)HH(iqkRx)-2mHm{iqkRi)}  ]  =  0 


Q 


(4.14) 


(4.15) 


From  eq.  (4.14)  and  eq.  (4.15),  one  dispersion  relation  is  derived  as 


1         Hm{iqK)       2m      J  \iqK)  2m 
a      nHm(iqK)    {*_m)nJm  {°lm) 

{(--mf-4}[UHm{iqK)-     2m    }-{J^iqK)-     2m    })  =  Q 
Q  llaXHm(iqK)    («       '    KJm{iqK)  'J 

Q  Q 


(4.16) 


where  a  =  —,K  =  kR  ,  q- 


(  m) 

Q 


R  is  the  jet  radius  at  the  interface,  J  and  H  are  general  Bessel  functions. 
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4.1.1  Axisymmetric  Disturbances  (m  =  0) 

Two  dimensional  disturbances  are  presented  in  r  and  z  coordinate.  The  dispersion 
relationship  from  eq.  (4. 1 6)  becomes  j 

•AOWw        H^K),  ,  1  H^iqK) 


a        { co 
Q 


-4. 


V0(i**)' 


=  0 


where  g  = 


2  -4 


The  asymptotic  form  of  dispersion  relationship  is  reduced  to 

(^2-4){l--^-  +  i(l  +  -i-)}  +  (l-l)^  =  0 
2qK    a      2qK  a 


(4.17) 


(4.18) 


4. 1 .2  Non-axisvmmetric  Disturbances  (m  g  0) 


The  dispersion  relationship  for  m*0is  eq.  (4.16).  There  are  two  kinds  of 
instabilities.  One  is  two-dimensional  disturbances  in  case  of  K  =  0 .  Another  is  three 
dimensional  disturbances  in  case  of  K*0.  For  K  =  0,  two-dimensional  disturbances 
(azimuthally  periodic  disturbances)  are  presented  in  r  and  9  coordinates.  For  K  *  0 , 
three  dimensional  disturbances  (helical  disturbances)  are  presented  in  r,  6  and  z 
coordinates.  Although  there  is  no  axial  velocity  in  two-wheel  flows,  these  flows  might  be 
considered  as  infinite  cylindrical  flows  for  the  disturbance  in  z  coordinate. 
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For  K  =  0 ,  the  asymptotic  form  of  the  eq.  (4.16)  becomes 


1 


2m    ,  .         2m  . 
-}  {m  +  ^r,  }  + 


(  \.){m- 

a  ,G>      x  ,a> 


{(--«)  -4}[-{w  +  —  }  +  {™-  —  

(  m)  (  m) 


(4.19) 


}]  =  0 


For      0 ,  the  dispersion  relationship  in  asymptotic  form  of  eq.  (4.16)  becomes 


{(—  -m)  -4}-[  (l-a)  +  \   {-  +  } 

O  (l  +  ay      '      +  J  X2  (*_m) 

Q 


I  +  (»_w)  (<*_m)2 

+  —  (32m4 -8m2  +18)}  ^-f— )  {I-i— 

Q 

{— (32/w  -8w  +  18)  +  T  +  —  +  —  }]  =  0 
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4       /®  N       ,®  x2 

(  m)    (  m) 


(4.20) 


4.1.3  The  Effects  of  Surface  Tension 


In  Reshotko  and  Monnin's  mathematical  theory,  there  is  no  consideration  of 
surface  tension.  This  work  extended  Reshotko  and  Monnin's  theory  to  include  the  effect 
of  surface  tension.  When  we  consider  the  effects  of  surface  tension,  the  dispersion 
relationship  for  m  =  0  becomes 


-[{\  +  T(l-K2)}-a){- 
a  ( a 


Hx{iqK)^  ,  1  Hx{iqK) 


J0(iqfO 


H0(iqK)'  aH0(iqK) 


JJjqK) 
J0{iqK) 


0 


The  asymptotic  forms  of  dispersion  relationship  for  eq.  (4.21)  is  reduced  to 
{co2-A){\ 


-^-  +  -(l  +  -^-)}  +  [-{l  +  r(l-^)}-lM  =  0 
2qK    a      2qK  a 


the  general  dispersion  equation  for  m  *  0  becomes 
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(4.21) 


(4.22) 


«»  -  w)»  -  4}[-{//"  ( W  -  -  Jg-|  -  { J-  <  W  -  =  Q 

UQ  tlaXHm(igK)    («       '    {Jm{iqK)  (^>J 

fi  Q 
For  the  asymptotic  forms  of  dispersion  relationship  for  m  *  0  and  AT  =  0  becomes 

[-{1  +  7(1  -  m-)}  -  l){m  -  -^—}{m  +  -^-}  + 

f,ffl  >1  .,,1,  2wj  -  .  2/w  „ 
{(---)  -4}[-{W  +  7F— :}  +  {^--  }]  =  0 


(4.23) 


(4.24) 


(— -m) 


(--w) 
Q 


For  the  dispersion  relationship  for  m  *  0  and  A;*0,  the  dispersion  relationship  in 
asymptotic  form  of  eq.  (  4.23)  becomes 
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,(_-»)  -4)-[.+  — 


(l-m2  -K2)]- 


_jK_ 
(l  +  a) 


(!-«)  + 


1-«V,1  2m 
{-  + 


/  vi  +  C/;0a:14 


2m  4/77" 

+  ■ 


1    fl-cn  ,1       2m  , 
b  +  —  } 


l  -a2 

+  (32m4  -8m3  +18)} 

128  (gK)2U  +  aJ  h  {<o_m) 

r  1  4     o    2     io\      1  2/W  4m3  A 

{ — (32m  -8m-  +18)  +  -  +  +  }1  =  0 

(  m)    (  my 

Q  Q 

In  Chapter  6,  the  special  case  of  m  =  0  will  be  presented. 


(4.25) 


4.2  Analytical  Solutions  for  Annular  Fuel  Sheet  Moving  Axially  and  Rotating 

Tangentially 

The  modified  Bessel  function  of  eq.  (4.9)  is  applied  here.  To  obtain  the  radial 
disturbance  amplitude  u  in  terms  of  modified  Bessel  function,  let 

c,  =ipjQj\4-nJ2)Ak-m   and  c2  =  i/>/i.2(4-  n/)Bkm  set.  From  Eq.  (4.9),  Dp  is 

obtained  as 


Dp  =  ciIm(kqJr)  +  c2Km(kqjr) 


(4.26) 


From  eq.  (4.5),  u  in  terms  of  Dp  and  p  is  also  expressed  as  the  combination  of  the 
first  and  second  kind  of  modified  Bessel  function.  Therefore,  u  is  expressed  as 

N. 


u = kqjNj  [Ak-»  ^1 + ^}     + Ens, K;    + — }  ^4^] 


or 
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where       =  kq}r 


(4.27) 


The  coefficients  of  A  and  B  in  eq.  (4.27)  are  defined  as  and  B0  on  the  axis. 
Similarly,  the  coefficients  of  A  and  B  are  defined  as  Al  and  5,  within  inner  air  layer,  A2 
and  B2  within  fuel  layer,  A3  and  53  within  the  outer  air  layer. 

u  is  finite  at  both  the  axis  and  the  outerwall.   When  £ .  is  close  to  the  axis, 

Km(€j)  8oes  t0  00  •  Thus>      should  be  0  and  Ax  will  be  unknown.  On  the  other  hand, 

LiCj)  goes  to  oo  as  £y  approaches  to  the  outerwall.  Therefore,  A3  should  be  0  and 

#3  will  be  unknown.  Ultimately,  there  are  4  unknowns  such  as  A\  A2  B2  and  B3  •  To 

find  the  4  unknowns,  four  equations  are  needed.  Four  equations  are  consisting  of 
interfacial  conditions  for  the  displacement  and  pressure  on  the  interface  on  both  sides  of 
annular  fuel  sheet.  Thus,  two  equations  are  described  at  inner  interface  of  fuel  sheet  and 
outer  interface  of  fuel  sheet,  respectively.  First,  eq.  (3.2.4)  and  eq.  (3.2.1 1)  are  applied  to 
the  interface  between  the  inner  air  and  fuel  sheet.  The  interfacial  conditions  become 


(4.28) 


and 
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4  (<?2^.)  «2 

t 

a^BT [W  -  4)  +  {q^K:  +  )  =  0 


(4.29) 


where  JC,  =  kRx ,  a2i  =  —  and  /?,,  =  — -2- 


When  the  surface  tension  is  considered,  eq.  (4.29)  becomes 

^r"[(V-4)+{^f|^  +  ^)]/.(^,)- 

/ 


«2,A,-^"m[(«2  -4)  +  f 
«2,A1252r[(/J22-4)  +  { 


A,  (?2*l)  «2 

q^Kj  (q,^)  ,  2m 


+  — }{l  +  7;,(l-^-^)}]/m(^^)- 


(?2*l  )  » 


+  —}{!  +  721  (1  -  m2  -  K,2 )}]     (q2Kt )  =  0 


where  7Z,  = 


AA 


(4.30) 


Second,  eq.  (3.2.4)  and  eq.  (3.2.1 1)  are  also  applied  to  the  interface  between  the  outer  air 
and  fuel  sheet.  The  interfacial  conditions  becomes 


2      L     I  (a  K  \     +—iJ^2^2)+B2k  [  + — ]Km(q2K2)- 

Jm^2K2)        n2  Km  (q7K?)  n7 


B3k™[ 


^K2Km  (q,K2)  ,  2m 


(4.31) 


Km(q3K2)  n 


+—]Km(q3K2)  =  0 


and 
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k«22  -  4) + &-K;Jr  {l\2) + — >i    ) + 

B2m»22  -  4) + {q'Kf:  {qv? + — >]f-(^^ )  - 


^m  ($3*2 )  "3 


where  K2  =  kR2 


(4.32) 


n,  = 


1  *"i 
N,     kW^  +  mQ,  -  co 


"3  = 


9,  = 


?2  = 


9:, 


4 

•4 

V"32  - 

4 

When  the  surface  tension  is  considered,  eq.  (4.32)  becomes 


A2k-[(n>  -  4)  +  {q>Y;  (^f2)  +  — }  {1  I  r22(l  -  m-  -  K22)}]Im(q2K,)  + 
52r  [(«22  -  4)  +  +  ^}  {1 1  T22(l  -  m-  -  K2>)}]Km(g2K2 )  - 

Km  (<?2A2)  «2 

Am  W3^2)  "3 


(4.33) 
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where  a„  =  — ,  /?„  =  —  and  7U  =  — r . 

From  eq.  (4.28)  through  eq.  (4.33),  there  are  four  unknowns.  The  determinant  should  be 
zero  to  solve  the  equations. 


I 22  K2 


K 


3: 


32 


0 

^2 

0 

0 

^4 

0 

where  Xt  =  A^",  X2  =  A2k-m,X3  =  B2km,  X,  =  £3£m  and 


(9:*:) 

Km(<llK2)  "3 


/u4=[(V-4)  +  {^4f^  +  ->W) 


(4.34) 


/2I*  =  a2  A2[(«22  -  4)  +  {^r-^  +  —  >]/„<*,*, ) 

L(^2K\)  n2 

V  -[(«.■  -4)+{24f^  +  ^)]/.(?,^) 
Im(q2K2)  nx 
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Km^2K\)  n2 


K2;  =[{n*-A)  +  { 


q2K2Km  (q2K2)  ,  2w 


(4.35) 


Km{q,K2)  n, 


when  the  surface  tension  is  considered,   I2*,I22',K2'  and  #22*  in  eq.  (4.35)  become 


/„*  =  a2i /?212[(«22  -  4)  +  (^f[)  +  ^} {1  +  rai (1  -      -  KfWJq^ ) 

I 

C  =  [(n22  -  4)  +  &-K;J;  (l\2)  +—}{!  +  ^(1  -  m1  -  K2>)})Im(q2K2) 

L^2K2)  n2 


K2;  =  aMUn*  -  4)  +  (&&£i&5>  +  — }  {1  +  rai(l  -  *r  -  K,2  JSC, ) 

KM2K\)  n2 


K22'  =  [(«22  -  4)  +  {?2^"  +—}{!  +  i(l  -  ™2  -  K2*)}]Km(q2K2) 

Km{q2K2)  n2 


(4.36) 


where       a-,,  = 


.A 


,a32  = 


'/^32  ~~     ^        '  ^2!  — 


2/ 


p2Q22/?,3 


and 


r  = 


p2q22r2' 


When  the  determinant  of  the  coefficient  matrix  is  set  to  zero,  eq.  (4.34)  becomes 


Al^22^2l  -^32    ~  Al^32^21  "^22    —  A  1^21  A>2  ^32    +  A  1  ^32  ^22  ^21  ~ 
4l  ^22  A  I  ^32    +  4 1^32  A  I  ^22    +  ^22^21  Al  ^32    ~~  A  I  ^32^22  ^21     =  0 


(4.37) 


/ 
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Here  the  first  derivatives  of  modified  Bessel  functions  are  expressed  as 


m 


L  (9,^)  =  /.+,  (qJKJ)  +  —-Im(qJKJ) 


(4.38) 


m 


(4.39) 


Km  (<7,K,)  =  -Km+1  (qJKJ)  +  —-Km{qjKJ) 

qjfi-j 

Eq.  (4.37)  is 


the  type  of  F(^,^,^,^,qi,q,,q3,ni,n,,n3,K]iK2,m,  —  ,—)  =  0 
p2  p3  Q,   fi,  ii2  S23 

where  co  is  the  only  unknown  parameter  to  be  solved.  Since  the  complex  eigen 
amplification  factor  falls  inside  the  arguments  of  the  Bessel  functions,  general  features 
of  the  system  are  difficult  to  see  without  the  use  of  a  computer.  However,  some 
important  results  can  be  obtained  by  using  an  analytical  expansion  for  sufficiently  large 
or  small  axial  wave  numbers.  In  those  cases,  the  Bessel  functions  reduce  to  an 
exponential  form  or  a  power  series. 

Since  eq.  (4.34)  consists  of  I  and  K'of  Bessel  functions,  I  and  K  can  be 
transformed  into  asymptotic  expansions. 

When  co  is  the  conjugate,  disturbances  are  growing  or  decaying  according  to 
whether  col  is  positive  or  negative,  respectively.  Since  the  fluids  have  been  assumed  to 
be  inviscid,  negative  mt  for  damping  the  disturbance  is  not  considered. 

The  physical  phenomenon  of  fuel  sheet  breakup  in  gas  turbine  combustor  is  due  to 
the  instability  of  annular  fuel  sheet  moving  with  axial  and  tangential  velocity  in  the 
surrounding  gas  which  has  axial  and  tangential  velocity.  Helical  flow  annular  jet  includes 
the  flow  characteristics  of  both  the  axial  flow  annular  jet  and  swirling  flow  annular  jet. 
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Before  investigating  the  physical  phenomena  of  the  annular  fuel  sheet  axially 
moving  and  tangentially  rotating  exposed  to  the  surrounding  air  axially  moving  and 
tangentially  rotating,  we  first  investigate  independently  two  kinds  of  instability  included 
in  the  physical  phenomena  of  the  annular  fuel  sheet.  First  kind  is  the  instability  for  the 
annular  sheet  moving  in  axial  direction  in  the  surrounding  gas  moving  in  axial  direction. 
Second  kind  is  the  instability  for  the  annular  sheet  rotating  in  tangential  direction  in  the 
surrounding  gas  rotating  in  tangential  direction. 

The  degrees  of  these  instabilities  are  basically  influenced  by  the  axisymmetric 
mode  (m  =  0  and  k*  0)  or  many  azimuthal  modes  (m*  Oand  k  =  0)  or  many  helical 
modes  (m*  0  and  k*  0). 

The  axisymmetric  mode  for  m  =  0  is  limited  to  two-dimensional  disturbances  in  r 
and  z  coordinates.  Therefore,  this  mode  is  independent  of  6 .  But  the  component  of 
perturbation  velocity  v  is  not  equal  to  zero.  The  azimuthally  periodic  mode  for  k  =  0  is 
limited  to  two-dimensional  disturbances  in  r  and  0.  Therefore,  this  mode  is 
independent  of  z  .  The  helical  disturbance  mode  for  m*  0  and  k*  0  is  three-dimensional 
disturbances  in  r  ,  0  and  z .  The  properties  of  the  instabilities  to  be  investigated  is  as 
following: 


Table  4.1  The  special  cases  for  which  analytical  solutions  will  be  derived. 


Axial  flow  annular  jet 
Section  4.2.1 

m  =  0  and 

k*0 

Rotating  flow  annular  jet 
Section  4.2.2 

m  =  0  and  k*  0 
Section  4.2.2.1 

m*  0  and  k  =  0 
Section  4.2.2.2 

m*  0  and  k*  0 
Section  4.2.2.3 

Helical  flow  annular  jet 
Section  4.2.3 

m  =  0andk*0 
Section  4.2.3.1 

m*  0  andk  =  0 
Section  4.2.3.2 

m*  0  and  k*  0 
Section  4.2.3.3 
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4.2.1  The  Axial  Flow  Annular  Jet 


There  is  no  rotational  velocity  component  in  flow  region.  Therefore,  Q,  Q2  and 
Q3  are  all  zeroes.  Also,  m,  circumferential  wave  number  is  not  considered.  Eq.  (4.35) 


becomes 


/2,=i 


/„  =  K}2=0 


I,  =  "I 


K, 


*o'(*,) 
/„"(*.) 


K2]  -  - 


K0(K2) 


I22  =  -p2N{ 


K22  =  -p2N2 


2K0{K,) 


*»  ~  2  u^y 

^33  =  PlK~ 


(4.40) 


K}:  =  p3N32 
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When  the  surface  tension  is  considered,    I22 ,  I22 ,K22 ,K22'  and  Ki2  *  in  eq.  (4.35) 


becomes 


U  =-p.N.2{^±r  +  K]  ^—(l-Kf)^^-  } 


i'  =  p.n:-  {i + tL  — ^-r(i  -  k;-)^^-  } 
-    -  ■      P2W     '  h(K2) 


22     ^2  2  l/0(*2)      ~p2N22R23"  2JI0(K2) 


3"        Pl     2   V2^22    /0  (*2>  *  P2^2^2  h  (K2) 


(4.41) 


The  asymptotic  form  ofeq.  (4.37)  becomes 


Ml 

P2^ 


- 1}  0  -  ^ttt}]  -  exp2^[{^  +  1}{1  + 


PA 


Ml 
P2K22 


}]  =  o 


(4.42) 


where 


N,  =kW7-co 


N3  ^kW^-co 
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K,  =  kR, 


£  = 


R2-R, 
R 


When  the  surface  tension  is  considered,  Equation  (4.42)  becomes 


l}{l-^V}]-exp2^[{^  + 


Mi 


(4.43) 


where 


Y,  =K.  V^(l-^,2) 


Equation  (4.42)  is  called  the  dispersion  relationship  for  annular  axial  jet. 
When  p]  =  p3  and  N{  =  N3 
eq.  (4.42)  becomes 


[(^4  - 1}  - 1}]  -  exp'-If^  -  1}  (1  +  =  0 


(4.44) 


and  eq.  (4.43)  for  the  effects  of  surface  tension  becomes 


■  1}]  -  exp2^[{};  -  Ml.  _!}{!  +  Mil}]  =  o 


p.TV,2  p2N{ 


pA\ 


P2*2 


(4.45) 
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4.2.2  Rotating  Flow  Annular  Jet 


There  is  no  axial  velocity  component  in  flow  region.  Therefore,  Wt ,  W2  and  W3 


are  all  zeroes.  Equation  (4.35)  is  applied.  Therefore,  each  term  in  the  coefficient  matrix 
in  eq.  (4.34)  becomes 


q2KJm  {q,Kx)  2m 


qiK2l :^  +  2m  ^ 


K2]=[ 


L  i^K2)  n2 
q1K1K1J32K1)2m 


Km(q2K,) 


+  —]Km(q2Kl) 


q2K2Km(q2K2)  +  2m 

Km(q2K2)        n2     fflW"  2) 

giK2Km' (giK2) +  2m 
Km(q,K2)  n/ 


Ai  —  LC«i  -4)  +  t — — — 77T~  +  — UymWiAi) 

/22*  =  [(V  -  4)  +  {^f^  +  ^}]Im(q2K2 ) 


k2;  =  «21/?2,W  -4)+{^§r  y,} +—}]*.(?,*,) 

k2: = -  4) + + ^k<*A) 

Km(l2K2)  n2 
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k2; 


*»PM  ~  4)  +  &-Kf"  ^  +  -}]^*, ) 

Km(<}2K\)  n2 

=  iin,;  -  4)  +  {^-'y^  +  M]*.<*A) 

Km(<l2K2)  n2 


(4.46) 


32 


= ^  [(>v  -  4) + {9>Yr  {V  f2) + -}]^ (^ } 


where  a,. 


fh 


'  ^"32  — 


.A 


32 


2J 


and 


"i  = 


n3  = 


?i  = 


q2  = 


Nx  _  wQ,  -  o) 

7V2  _  ffjQ,  -  CO 

Q2  "  Q2 
7V3  _  wQ3  -  <y 
Q3  "  Q3 


vwi2  - 

4 

V«22  ~ 

•4 

"2 

4 

when  the  surface  tension  is  considered,  721  ,/22  ,K2l  and  £22  in  eq.  (4.46)  becomes 


A.^a^r^-^  +  ^'^^f'^—  }{l  +  r2l(l-^-^r)}]/m(^,) 

t 

I2;  =  [{n;  -  4)  +  (ft*,)  +       {1  +  r  0  _  m2  _  K;)}VM2K2) 

4(^2^2)  W2 
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K2;  =  a„ /?2,2[(«22  -  4)  +  {^f"  &ftf  +  — } {1  +  Tn (l-m2-  tf)>K.<fc*i ) 


Kv'  =  [(«22  -  4)  +  {^Mm  (^2)  +       {1  +  r  0  _  m2  _  ^^(g a^ ) 


where 


(  \ 


a- 


\P\) 


•fin 


,a 


32 


A. 


r  = 


and 


(4.47) 


T  = 


p2Q22&3 


4.2.2.1  Axisvmmetric  disturbances(m  =  0  and  k  g  0) 


When  eq.  (4.37)  is  expressed  in  asymptotic  form,  eq.  (4.37)  becomes 

-*M**Pnffi% exp-^(K22  -4  +  L2^)(»32  -4-^2) 
-«2. A,2«32 A:' exp^(«22  -  4  -  J  tf,  )(V  -  4  -  ?3tf2 ) 


+«2,^,2)/%^%  exp-^(«22  -4  +  ^2/:,)(«22  -4-^2) 
-«2,A,2 exP?2V(^2  -4  +  ^2)("22  -4-^,) 
+«32^22^j%  exp-^(n,2  -4  +  i*j XV  -4-^3^) 

-«32/?322  ^^^^  «P«*V  -4  +  ^kX"32  -4-^2) 

-^exp-<^(n22  -  4  -  <72£2)(«,2  -  4  +L,  Aj) 


?3> 


exp^>22  - 4  +  <72tf2)«  - 4  +  4,*,)  =  0 


(4.48) 


When  the  surface  tension  is  considered,  the  asymptotic  form  of  eq.  (4.37)  becomes 
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[n22  -  4  +  q2K, {1  +  r„(l  -  K2)}](n2  -  4  -  93K2) 
[n22  -  4  - {1  +  r21  (1  -  K;)}]{n,2  - 4 - tf3£2) 

[«22  -  4 + ?2#,  {i  +  r2l(i  -  K2)}}[n2  -  4  -  q2K2  {i  +  r22(i  -  k22)}} 
[«22-4  +  <72/c2{i  +  r22(i-/:22)}][Aj22-4-e2^{i  +  7;,(i-/:,2)}] 

«p**V  ~ 4 - ft***1  +  rao - ^2)}](«,2 - 4 + ) 

-  j%  exp^[«22  -  4  +  42K2{1  +  7/22(l  -  X22)}](«,2  -  4  +  q,Kx)  =  0 


4.2.2.2  Azimuthallv  periodic  disturbances  (m  *  0  and  k  =  0) 


(4.49) 


The  asymptotic  form  of  eq.  (4.37)  becomes 


nx(nx  -2)-ar21/7,,2M,(tt2  -  2)  +  m(l  -  a21/?212) 


«,(«,  -2)  +  a21y?212«2(n2  +2)  +  m(l-a2]^212)V^2y 

h2(h,  -  2)  +  <?32/?322ft3(ft3  +  2)  +  m(l  -  avJi2)  _ 
n2{n2  +2)-a32pi22n3{ni+2)-m{\-av_pi2) 


(4.50) 


When  the  surface  tension  is  considered,  the  asymptotic  form  of  eq.  (4.37)  becomes 
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a*P*"*Pii2 0"  +  — X-w  +  — )K2  -  4  +  (m  +  —  ){1.4T21  (1  -  m2 )}] 


(«3-  -  4  -  m  +  — ) 


UJ 


2m_  2m 


2m, 


«2i  Ai  «32 (OT  +  — )(«  +  — )K  -  4  +  (-iw  +  — ){\.+T2]  (1  -  m" )}] 


n, 


n-, 


ru 


,  2  a  2m. 
(«3_  -  4  -  m  +  — ) 


2m 


2m 


2m, 


-a21/?2,2  (m  +  — )(-m  +  — )[n22  -  4  +  (m  +  — ){l+T2] (1  - m2)}] 


ru 


[n22_4+(_OT+ — ){l.+7i2(l-m2)}] 

»2 


/       \  2m 


2m 


2m 


2m, 


tf21/?212(m  +  — )(-m  +  —  )[>v  - 4  +  (m  +  — ){1.+7;2(1  -  m')}] 


ru 


2m 


[^_4  +  (_w  +  _){L+7;i(l_^)}] 


.  2.  2mw  2mw  .  2wr 
-a32J3i2-(m  +  — )(-m  +  — )(«,"  -  4  +  m  +  — ) 


Hi 


/  2  ,  2m 
(«3  -  4  -  m  +  — ) 


2m 


2m, 


2m, 


ai2fii2-(m  +  — )(-m  +  — -  4  +  m  +  — ) 


n. 


,       2         „  2WN 

(/?3_  -4-m  +  — ) 


2mw  2m. .  ,  .  2/rr 
(m  +  — )(-m  +  )(«,"  -  4  +  m  +  ) 


Hi 


2m, 


[«,2  -  4  +  (-m  +  ){l.+r„(l  -  m2)}] 


2m 


2m 


'V2 

2m, 


(-m  +  )(-m  +  — )(«,  -  4  +  m  +  ) 


2m 


4  +  (m +  —){!.+ r22(l  -  m2 )}] 


(4.51) 


4.2.2.3  Helical  disturbances  (m*  0  and  0) 


The  asymptotic  form  ofeq.  (4.37)  becomes 


a2Xfi2x2anpi22  exV^K*\qxKx  +  ~)(-q2K2  +  — )(«32  -  4-ftJC,  +— ) 


[n22-4  +  (q2K]+—)) 


[n2'-4  +  (-q2Kl+—)] 

-a2l/32>  WK^\qA  +^)(-q3K2+^)[n2>-4  +  (q^  +— ) 
[n22-4  +  (-q2K2+—)] 
a2x/32?  exp*^'^,*,  +—X-q3K2  +  l%«22  -4  +  foaJC.  +  — ) 
][»,»-4  +  <-gaAj+2*)] 


f   ->    a       v  2m 
"3 

2  »v„(?:*:-»2*i) 


2m 


2m, 


2m, 


+a32fi32-  cxp^^K'\q2K2  +  —)(-q2K]  + — )(«,  -4  +  q]Kl  +—) 


fU 


/  ■>  a  is  2m. 
(n{-4-q3K2  +  ) 

"3 

(o,/c,-o,a:,)/    v     2m              2m      -,  2m . 

+  expw-  1  92  -'fa,  AT,  +  )(-q3K2  +- — )(«,"      +        +  ) 


2m. 


[«22-4  +  (-^2+— )] 


2m, 


[/v-4  +  fo2*2  + — )]  =  0 


(4.52) 


When  the  surface  tension  is  considered,  the  asymptotic  form  of  eq.  (4.37)  becomes 


*ixPj**fia  exp^'-^fo,*,  +—  )(-q2K2  +  ^)(n2-4-q3K2  +—) 


n. 


2m 


[n22  -  4  +  (q2K,  +  — ){1  +  T21(l  -  m2  -  K2)}] 


-a^JanPn  +  — X<72*2  +  — )(»32  -4 -?,tf2  +  — ) 


n, 


-  4  +  (-^  +  —  ){i  +  r2l  (i  -  m2  -  k,2 )}] 

-a2l/?212  exp^^^K,  +  ^)(-^2  +  —){n2  -  4  +        +  — ) 


2»r 


{l  +  r21(i -  ™2  - ^,2)}[»22  - 4 + (-«72*:2  +  — ){i  +  r22(i - m2  - k2)}} 


a2,/y  exp^-^.K,  +  —)(-q3K2  +^)[n2-AHq2K2  3) 


2m, 


{1  +  r22(l  - m2  -  K22 )}][«22  - 4  +  {-q2Kx  +  — ){1  +  7/21  (1  - m2  -  K2)}] 


-ap/V  exp^^K,  +— X"<72*2  +  —)(n2-4  +  q^  -3) 


y   f     m       j,  2m 
(ni--4-qiK2+  ) 


+a320322  cxP^--^\q2K2  +  — X-^  I  — )(«,2  -  4  +  qxKx  +  ^L) 


/  ■>    ^       w  2/n 

(A73--4-^3^2+  ) 


+  exp 


2m, 


[n2  -4  +  (-q,K,  +  ){1  +        -  m2  -  K2)}} 


(4.53) 


-exp 


2  q^\-q2K,  +  X-ft*2  +  KK  -4  +  9,^,  +  ) 


2m, 


[«22  -  4  +  (q,K,  +  ){1  +  T„(l  -  m2  -  K2)}}  =  0 
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4.2.3  The  Helical  Flow  Annular  Jet 


4.2.3.1  axisymmetric  disturbances  (m  =  0  and  0) 

m,   circumferential  wave  number  is  not  considered.   Each  term  in  eq.  (4.34) 


becomes 


4)  (^2*2) 


K2I = [ 


Ao(«2*|) 


q2K2K0  (q2K2) 

A22  =  L  „  ,      „  .  J*oW2A2) 

32  ~  1  Y  In  Y  \  1KM>K2> 


/„*=[(V-4)  +  {^f^}j/„(^) 

721-  =  a2,/?212[(«22 -4)  +  {^^o(^fl)}]/0(^,) 

C=[(V-4)  +  {^/o(^f2)}]/0(^2) 
4(02*2) 

=  «2,  a,2[(«22  -  4) + iq"-Kf;  (*f 'Wft*.  ) 

^0(^2^1) 

iCn*  =  [(«22  -  4)  +  {M^MMl2}]ic  (^^} 

^0(^2^2) 


(4.54) 
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where 


n, 


N]  _kW]-co 
Q,  "  Q, 
N2  _kW2-a) 
Q2  "  Q2 

N3  _kWi-CO 

F^ 


12  = 


03 


fF 

fF 


(4.55) 


When  the  surface  tension  is  considered,     ,/2      and  /C,2  in  eq.  (4.55)  are  replaced  by 

n2KX(q2KA 


V  =  «2A2[(«22-4)+{ 

4*=[(«22-4)  +  { 


?2k2/0  (<?2a:2) 


/0(*2*2) 


}{i  +  r22(i-^)}]/0(^) 


,  ^ =[(^ -4)+{^°  y^^o+r^a-^2)}]^^)  (4-56) 

K2;  =  a,J2;[{n2  -  4)  +  {*2^*°  (*2f f}{l  +  r2I(l  -  K;)))KQ{q2Kx) 
When  eq.  (4.55)  is  changed  into  asymptotic  expansions  ,  eq.  (4.55)  becomes 


-«2.A.2«32^322  )/^^^5^eXP 

[n22-4  +  q2Kt](n2-4-q3K2) 


q2Kxf-  '/qyK2**P 


-«2,/?2,2«32/?32Y *Y  -  ^ 

[n22  -4-q2Ki](n32  -4-^2) 
+«2.P2i  y  y  exp 
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[n2  -4  +  q2K^\[n2  -4-q2K2] 


n  R  2  fatK\/  m^i7  Pvn(^2 
[n22  -  4  +  q2K2][n22  -  4  -  q2Kt] 


(n2  -4  +  q}  K,){n2  -4-q}  K2) 
(n2  -4  +  qiK])(n2  -4-q3K2) 


[n22 -4-q2K2](n2 -4  +  qxKt) 

/q^yj  /q2K2™> 
[n22-4  +  q2K2](n2-4  +  qlKi)  =  0 


When  the  surface  tension  is  considered,  eq.  (4.66)  becomes 


[n2  -4  +  q2Kx{\  +  7/21(l  -  K2)}]{n2  -  4  -  q3K2) 
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n    ft  2  \^-K\/        klKl/  exn 

(«,2  -A  +  q,Kx){n2  -4-q3K2) 

v  AM  A*K *  v 

[n22  -  4  -  q2K2  {1  +  T22(\  -  K22)}](n2  -  4  +  qxKx ) 

*^Jg3K/q2K2 
[n22  -  4  +  q2K2 {1  +  T22(\  -  K2  )}](n2  -  4  +  qxKx )  =  0 


(4.58) 


For  non-axisymmetric  disturbances  ( 0),  eq.  (4.35)  becomes 
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K2;  =  aa [(«22  -  4)  +  {*K*K:  TV1  +  > 

i 

when  the  surface  tension  is  considered,  //J-^/^'and  A,,*  in  eq.  (4.59)  are  replaced  by 
V  =  [(V  -  4)  +  (52f'}  +—}{!  +  r2, (1  -  m2  -  K,a))]/.(ftAj ) 


V  =  [(>V  "  4)  +  (^2f2)  +—}{!  +  r22(l  -  m-  -  K22)}]Im(q2K2 ) 
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4.2.3.2  Azimuthally  periodic  perturbations  (m  *  0  and  k  -  0) 


When  eq.  (4.37)  is  expressed  in  asymptotic  expansion  forms,  eq.  (4.37)  becomes 

_  ,      n  ->/       2mw         2m. r   2         /  2m 
a2ifi2\'arJn~(m  +  — )(~m  +  — )K"  -  4  +  (m  +  — )] 


/    1  2m 

(«3  -4- m  + — ) 


u2. 

_  2mw      2mxr   2        /  2m 

-ai\Pi{ar.Pv:{m + — )(w + — -  4 + + — )] 


/    1      ^  2m 

(«3-  -  4  -  m  +  — ) 

_  2/      2m.,        2mNr   2    >    /  2m 
-a„^„- (m  +  — )(-in  +  — -  4  +  (m  +  — )] 


r    2     a     /  2m 
|V  -4 +  (-*?!  + — )] 

«2 


/       \  2m 


_  2/  2m.,  2m. .  ,  .  ,  2mvl 
a2l02l2(m  + — X-W  +  — )[n2-  -4  +  (m  +  — )] 


r   ■>    a    /  2m 
[<-4  +  (-m  + — )] 


_  2.  2mN,  2mw  2  .  2m 
-a32j032  (m  +  — )(-m  +  — )«  -  4  +  m  +  — ) 


2m 


(n{  -  4-m  +  ) 


»3  v&J 

_  2/  2m..  2m..  2  2m. 
avpv-{m  +  — )(-m  +  — )(/?,-  -  4  +  m  +  — ) 

/  2  2m. 
(«3"  -  4  -  m  +  ) 

r  ->  a  <  2mN1 
[«2"-4  +  (-m  +  )] 


UJ 


2m          2m.,  2     .  2m. 
-(-m  +  )(-m  +  )(n,"  -  4  +  m  +  ) 


(4.61) 


n. 


0  yyi 

[n2 -4  +  (m  +  —  )]  =  0 


When  the  surface  tension  is  considered,  eq.  (4.61)  becomes 
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(4.62) 
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[n2  -  4  +  (m  +  — ){l.+r„(l  -  m2)}]  =  0 


4.2.3.3  Helical  disturbances  (m  g  0  and  k  *  0) 


Equation  (4.37)  expressed  in  asymptotic  form  is 
a,xpjanpn>  exp^'^,*,  +—)(-g2K2  +^)(n2-4-q3K2  3) 
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n2  n2  «, 
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"2 

When  the  surface  tension  is  considered,  eq.  (4.37)  becomes 
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(4.64) 


CHAPTER  5 

THE  VALIDATION  OF  NUMERICAL  FORMULATION 

The  validity  of  the  previous  numerical  formulations  is  investigated  in  this 
chapter.  There  are  three  known  analytical  solutions  available  to  be  compared.  Rayleigh 
(1879)'s  mathematical  theory  for  stationary  flow  assumes  an  inviscid  fluid  core  with 
surface  tension  surrounded  by  a  light  inviscid  fluid.  Reshotko  and  Monnin  (1965)'s 
mathematical  theory  for  rotating  flow  assumes  a  heavy  inviscid  fluid  core  with  rotation 
surrounded  by  a  light  inviscid  fluid  with  rotation.  Fung  (1983)'s  mathematical  theory  for 
rotating  flow  assumes  a  heavy  inviscid  annular  sheet  with  rotation  exposed  to  a  light 
inviscid  fluid  with  rotation.  The  computer  code  was  tested  against  these  known 
solutions.  In  addition,  Fung  (1983)'s  mathematical  theory  for  the  azimuthally  periodic 
mode  (m*  0  and  k  =  0)  with  no  surface  tension  was  compared  with  the  derivation  of  this 
dissertation,  eq.  (4.37).  After  the  computer  code  is  validated,  it  will  then  be  extended  to 
the  analysis  of  the  inviscid,  rotational  and  axial  flow  with  or  without  surface  tension. 

5.1  Rayleigh's  Stationary  Jet  with  Surface  Tension 

Rayleigh  (1879)  made  a  mathematical  analysis  of  the  instability  of  a  non-viscous 
jet.  He  considered  a  laminar  liquid  flow  with  a  velocity  potential  where  the  jet  was  only 
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under  the  influence  of  surface  tension  forces.  Although  the  conditions  considered  by 
Rayleigh  do  not  exist  in  an  actual  jet,  the  general  conclusions  of  Rayleigh's  work 
regarding  jet  instability  have  been  accepted.  In  his  analysis,  the  surface  waves  grow  as 


e°* ,  where 


r3p  i0m 


[l-rin 


(5.1) 


a  is  the  surface  tension,  p  is  the  density  of  liquid,  R  is  the  undisturbed  radius  of  the 
jet,  k  is  the  wave-number  and  I0 ,  7,  are  the  modified  Bessel  functions.  The  dispersion 

In 

curve  of  co  versus  k  shows  a  maximum  growth  rate  at  kR  -  0.697  or  X  =  ——  =  9.02  R . 

k 


It  is  verified  if  the  computer  code  shows  the  maximum  growth  rate  at  k  = 


0.697 
R 


This 


verification  is  shown  in  figure  7.  In  figure  7,  the  results  from  the  computer  code  shows 

0.697 


that  the  maximum  growth  rate  occurs  at  k  = 


R 


5.2  Reshotko  and  Monnin's  Two- Wheel  Flows  with  No  Surface  Tension 


Reshotko  and  Monnin  (1965)  investigated  the  instability  for  the  two-fluid  wheel 
flow  of  a  heavy  density  core  surrounded  by  a  lighter  fluid.  The  geometry  of  such  a  flow 
is  shown  in  figure  6.  The  flow  is  considered  to  be  only  azimuthal  and  solely  a  function  of 
the  radial  coordinate.  In  Reshotko  and  Monnin's  two  wheel  flows,  the  instability  of  two 
layers  with  the  same  rotational  velocity  and  different  density  was  investigated.  The  inner 
fluid  density  is  greater  than  the  outer  fluid  density.  For  a  heavy  fluid  core  surrounded  by 
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Figure  7.  The  effects  of  surface  tension  on  Rayleighjet 
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a  light  fluid,  the  inertial  nature  of  the  instability  is  apparent  from  the  pressure  distribution 
for  px>  p2.  The  heavy  fluid  encounters  a  pressure  gradient  inadequate  for  radial 
balance  on  penetrating  the  interface  and  the  element  of  heavy  fluid  will  continue  its 
motion  outward.  This  behavior  is  characteristic  of  Rayleigh-Taylor  instabilities. 

5.2.1  Axisymmetric  Disturbances  (m  =  0) 

The  computed  results  using  eq.  (4.17)  and  eq.  (4.18)  are  compared  with  the  results 
presented  in  figure  8  of  Reshotko  and  Monnin  (1965)'s  technical  note.  For  m  =  0,  the 
growth  rate  seems  to  quadradically  increase  as  axial  wave  number  K  increases.  It  is 
shown  that  the  computational  results  match  Reshotko  and  Monnin's  results  very  well. 

5.2.2  Non-axisymmetric  disturbances  (m  ^  0)  . 

Reshotko  and  Monnin's  results  in  figure  9  are  compared  with  this  works  solutions 
to  eqs.  (4.24)  to  (4.25)  in  figure  10.  In  figure  9,  the  growth  rate  at  m  =  1  seems  to 
increase  linearly  as  the  axial  wave  number  K  increases.  However,  the  growth  rate  for  m 
=2  to  m  =  5  seems  to  be  independent  of  K  when  K  is  equal  to  1  or  less  than  1 .  The 
growth  rate  increases  very  slowly  as  m  increases  and  K  increases  beyond  K  =  1 .  At  large 
K,  growth  rate  seems  to  converge  into  a  specific  value  independent  of  m. 

In  figure  10,  the  computational  results  are  shown.  The  growth  rate  at  m  =1  and  K 
=  0  is  slightly  below  Reshotko  and  Monnin's  result  while  the  growth  rate  for  m  =  2  to  5 


82 


The  comparison  between  numerical  code  and 
Reshotko's  results  for  axially  symmetric 
disturbances  for  two  wheel  flow 
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Figure  8.  The  growth  rates  on  m 


10 


12 


83 


case  :  two  wheel  flow 
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Figure  9.  Reshotko's  results  on  m  variations 


Figure  10.  computational  results  on  m  variations 
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and  K  =  0  is  similar  to  Reshotko  and  Monnin's  results.  As  K  increases,  the  growth  rate 
decreases  until  K  approaches  1 .  However,  the  growth  rate  increases  as  K  increases  over 
1.  There  are  some  differences  between  the  computation  code  and  Reshotko  and 
Monnin's  results  for  K  =  0  to  K  =  1.  But,  there  are  similar  trends  between  the 
computation  code  and  Reshotko  and  Monnin's  results  beyond  K  =  1.  Theoretically,  the 
growth  rate  is  expected  to  increase  as  K  increases. 

For  K  =  0  to  K  =  1  in  Reshotko  and  Monnin's  results,  the  growth  rate  does  not 
seem  to  be  correct  because  the  growth  rate  seems  to  be  constant.  On  the  other  hand,  the 
growth  rate  shown  in  figure  10  does  not  seem  to  be  correct  in  the  theoretical  aspect 
because  the  growth  rate  decreases  in  that  same  region.  In  private  communications  with 
both  Dr.  Reshotko  and  Dr.  Monnin,  co-authors  of  Reshotko  and  Monnin's  results, 
neither  had  a  clear  statement  of  how  different  results  were  obtained  when  their  equations 
and  those  used  in  this  dissertation  are  the  same.  The  computer  code  of  the  series 
expression  of  Bessel  fuction  was  testified  to  verify  the  validation  of  the  asymptotic 
approach.  However,  the  computational  results  are  not  produced  enough  to  validate  the 
asymptotic  approach.  Despite  these  differences  at  small  K  between  Reshotko  and 
Monnin's  work  and  the  present  work,  it  will  be  assumed  that  the  computer  code  is  valid. 

The  extension  of  Reshotko  and  Monnin's  work  to  include  surface  tension  was 
derived  in  section  4.1.3  of  this  dissertation.  The  prediction  of  the  effects  of  surface 
tension  in  Reshotko  and  Monnin's  equations  will  be  discussed  in  Chapter  6. 
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5.3  Fung's  Equation  for  Azimuthally  Periodic  Perturbations  On  *  0  and  k  =  0) 

with  No  Surface  Tension 


Fung  (1983)  established  the  general  dispersion  equation  for  three  rotating  layers 
of  fluids  as  following: 


-2)-a2lfi22n2(n2  -2)  +  m(\-  a2xP22) 
-2)  +  a2Xp22n2{n2  +  2)  +  m(\  -  a2xf322) 


2m 


n2(n2  -2)  +  ai2j3i22ni(ni+2)  +  m(l-ai2/3J22) 
n2(n2  +  2)  -  a32p32\(n3  +  2)  -  m(l  -  ai2Pn2) 


(5.2) 


CO 


where  n,  =  m-  —  ,  au  =  —  and  fiu  = 


This  can  be  compared  with  eq.  (4.50)  in  section  4.2.2  which  is  derived  with 
swirling  flow.  Equation  (4.50)  is  exactly  the  same  as  eq.  (5.2)  which  Fung  (1983) 
obtained.  Equation  (4.50)  is  one  of  several  solutions  to  be  derived  in  this  dissertation 
and  is  obtained  by  omitting  both  axial  annular  jet  and  surface  tension.  This  validates  the 
derivation  of  the  equations  used  in  current  computer  code.  Therefore,  it  is  assumed  that  it 
will  valid  when  the  current  code  will  be  extended  to  investigate  the  instability  of  the 


helical  flow  annular  jet. 


CHAPTER  6 
RESULTS  AND  DISCUSSION 

Although  the  surface  tension  was  not  included  in  Reshotko  and  Monnin's 
mathematical  approach,  the  effect  of  surface  tension  has  a  significant  effect  on  the  degree 
of  instability  between  interfaces  of  two  different  fluids.  The  computational  results  for  the 
effect  of  surface  tension  will  be  discussed. 

The  main  objective  in  this  dissertation  is  to  investigate  the  instability  of  an 
annular  fuel  sheet  with  axial  and  tangential  velocity  surrounded  by  the  air  with  axial  and 
tangential  velocity. 

The  computed  results  for  4  cases  will  be  discussed  in  the  following  order: 

1.  The  effects  of  the  surface  tension  using  the  extension  of  Reshotko  and  Monnin's 
equations 

2.  Annular  fuel  sheet  with  axial  velocity  and  surrounding  gas  with  axial  velocity 

3.  Annular  fuel  sheet  with  tangential  velocity  and  surrounding  gas  with  tangential 
velocity 

4  The  annular  sheet  with  axial  and  tangential  velocity  moving  through 
the  surrounding  gas  with  axial  and  tangential  velocity 
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6.1  The  Effects  of  the  Surface  Tension  Using  the  Extension  of  Reshotko  and  Monnin's 

Equations 

The  inner  denser  layer  and  outer  layer  rotate  in  a  clockwise  direction  at  the  same 
speed.  Equation  (4.22)  is  being  solved.  Figure  11  shows  that  the  growth  rate  first 
increases  and  later  decreases  with  increasing  axial  wave  number  K.  It  is  predicted  as  the 
surface  tension  in  the  inner  layer  gets  higher,  the  growth  rate  increases  at  very  smaller  K 
and  then  decreases.  This  predicts  that  the  disintegration  of  the  inner  layer  fluid  by 
aerodynamic  forces  is  more  difficult  with  higher  surface  tension.  This  can  be  explained 
from  the  relationship  between  the  pressure  difference,  surface  tension  and  curvature 
radius  at  the  interface  of  different  phases  in  eq.  (3.2.20).  The  wave  length  is  the  inverse 
of  the  wave  number.  The  radius  of  curvature  is  large  when  the  wave  length  is  large. 
When  the  curvature  radius  is  large,  the  effect  of  the  surface  tension  is  relatively  small. 
Therefore  pressure  difference  predominantly  depends  on  the  centrifugal  force  of  the  inner 
core  region  is  higher  than  that  of  the  outer  region.  The  instability  grows  at  the  interface. 
That  phenomena  leads  to  the  increase  of  growth  rate. 

However,  the  effect  of  the  surface  tension  is  relatively  large  when  the  axial  wave 
number  is  large  and  the  curvature  radius  is  small.  Pressure  distribution  depends  on  the 
difference  between  the  centrifugal  force  and  surface  tension  force;  therefore  the  degree  of 
instability  is  balanced  between  the  centrifugal  force  and  surface  tension  force.  The 
growth  rates  decrease  and  the  disintegration  of  the  inner  layer  fluid  gets  more  difficult. 
Thus,  the  growth  rate  increases  around  the  small  wave  number,  later  it  decreases  around 
the  large  wave  number.  Maximum  growth  rate  is  formed  at  a  specific  axial  wave  number 
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as  shown  in  figure  1 1 .  When  the  surface  tension  increases,  the  maximum  growth  is 
formed  at  smaller  axial  wave  number  since  the  growth  rate  increases  until  smaller  axial 
number  and  decreases  right  after  as  axial  wave  number  increases. 

6.2  Axial  Flow  Annular  Jet 

The  annular  sheet  was  injected  in  the  axial  direction  at  a  speed  different  from  the 
surrounding  gas.  The  principal  cause  of  instability  is  due  to  the  relative  velocity  between 
the  sheet  and  the  surrounding  gas.  The  sheet  interacts  with  the  surrounding  gas  whereby 
rapidly  growing  waves  are  imposed  on  the  sheet  just  after  the  nozzle  from  the  atomizer 
exit  of  the  gas  turbine  combustor.  The  velocity  jump  at  both  inner  interface  and  outer 
interface  of  the  annular  fuel  sheet  for  step-wise  approximations  of  continuous  functions 
yields  the  Kelvin-Helmholtz  instabilities  for  even  small  wave  lengths.  This  instability  is 
due  to  the  strong  shearing  rates  existing  in  the  sharp  shear  layer  created  at  the  vortex 
sheet.  Instabilities  produced  by  the  Kelvin-Helmholtz  effects  due  to  the  relative  motion 
at  interfaces  are  very  strong  and  cause  the  fuel  sheet  to  be  severely  unstable. 
Disintegration  occurs  when  the  wave  amplitude  reaches  a  critical  value  and  fragments  of 
annular  fuel  sheet  are  torn  off.  The  fragments  rapidly  contract  into  drops  under  the  action 
of  surface  tension.  Equation  (4.42)  is  solved  for  no  surface  tension,  eq.  (4.43)  is  solved 
for  the  effects  of  surface  tension. 

In  figure  12,  the  effect  of  surface  tension  in  the  annular  axial  jet  is  shown. 
There  is  a  significant  effect  of  reducing  growth  rate  at  high  values  of  surface  tension. 
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Therefore,  the  growth  rate  first  increases  and  decreases  later  as  axial  wave  number  K 
increases  and  involves  much  larger  surface  tension.  When  the  typical  value  of  surface 
tension  is  applied  to  the  analysis  of  instability,  the  growth  rate  decreases  little  from  the 
zero  viscosity  case. 

In  figure  13,  the  effect  of  fuel  density  change  for  the  annular  axial  jet  is  shown.  It 
is  shown  that  the  growth  rate  increases  as  the  axial  wave  number  increases  at  the  typical 
value  of  surface  tension  and  the  growth  rate  decreases  as  the  density  of  the  annular  axial 
jet  increases.  It  tells  that  the  mechanism  of  annular  denser  sheet  disintegration  should  be 
affected  to  some  extent  by  a  change  in  annular  fuel  sheet  density.  It  is  predicted  that  it  is 
harder  for  denser  annular  sheet  moving  in  the  axial  direction  to  be  broken  up.  For  the 
annular  fuel  sheet  of  high  density,  more  resistance  and  inertial  effect  to  disintegration  at 
the  atomizing  edge  is  expected  to  occur. 

In  figure  14,  the  effects  of  thickness  variation  for  the  annular  axial  jet  is  shown. 
It  is  shown  that  the  growth  rate  decreases  as  the  thickness  of  annular  fuel  sheet  increases 
but  the  growth  rate  did  not  decrease  any  more  beyond  a  certain  thickness  of  the  annular 
fuel  sheet.  This  similar  phenomena  is  described  in  J.  Meyer  and  D.  Weihs  (1987).  This 
result  shows  that  there  is  a  critical  penetration  thickness.  It  seems  that  the  annular  fuel  jet 
behaves  like  a  full  liquid  jet  when  the  annular  thickness  is  greater  than  the  critical 
thickness  of  fuel  sheet  and  the  growth  rate  is  independent  of  thickness.  When  the  annulus 
thickness  is  less  than  the  critical  thickness,  the  jet  behaves  like  a  two-dimensional  liquid 
sheet  and  the  growth  rate  increases  as  the  jet  thickness  decreases.  Therefore,  an  annular 
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Figure  14.  The  effects  of  thickness  variations 


fuel  sheet  with  a  sufficiently  small  ring  thickness  is  expected  to  disintegrate  into  drops 
much  faster  than  a  thicker  fuel  jet  disintegrates  into  drops. 

In  figure  15,  the  effects  of  the  axial  velocity  variation  of  the  annular  fuel  sheet  is 
shown.  It  is  shown  that  the  degree  of  instability  is  increased  as  the  axial  velocity  of  the 
annular  fuel  sheet  is  increased.  That  is,  the  growth  rate  is  increased  as  the  axial  velocity 
of  the  annular  fuel  sheet  is  increased.  It  is  expected  that  the  higher  the  relative  axial 
velocity  between  the  annular  fuel  sheet  and  the  surrounding  gas  is,  the  faster  the 
aerodynamic  forces  cause  the  annular  fuel  sheet  break  into  spherical  shells  with 
exponentially  increasing  amplitude. 

In  summary,  faster  breakup  of  the  annular  fuel  sheet  moving  in  the  axial 
direction  might  be  implemented  by  the  decrease  of  density  or  decrease  of  thickness  or 
increase  of  the  axial  velocity  of  the  annular  fuel  sheet. 

6.3  Rotating  Flow  Annular  Jet 

The  flow  which  is  now  considered  has  centrifugal  instability.  This  is  analogous  to 

the  Rayleigh-Taylor  instability  under  which  Reshotko  and  Monnin's  two  wheel  flow  is 
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classified.  There  is  the  annular  fuel  sheet  exposed  to  the  surrounding  gas.  The  annular 
sheet  is  rotated  with  the  angular  velocity  different  from  that  of  the  surrounding  gas.  The 
flow  regions  consist  of  three  layers  which  include  an  inner  gas  cylinder,  an  annular  fuel 
sheet  and  an  outer  annular  surrounding  gas  layer.  However,  the  angular  velocity  in  the 
outer  annular  gas  layer  is  considered  as  almost  zero  because  the  distance  from  the  center 
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of  the  combustor  to  the  outer  layer  radius  and  there  is  no  angular  velocity  in  most  of  the 
outer  layer  region.  Eq.  (4.48)  for  m  =  0  and  K  *  0,  eq.  (4.50)  for  m  *  0  and  K  =  0,  eq. 
(4.52)  for  m  *  0  and  K  *  0  are  solved  for  no  surface  tension.  Eq.  (4.49)  for  m  =  0  and 
K  *  0,  eq.  (4.5 1)  for  m  *  0  and  K  =  0,  eq.  (4.53)  for  m  *  0  and  K  *  0  are  solved  for  the 
effects  of  surface  tension.  Figure  16  and  17  show  that  surface  tension  does  not  affect  the 
growth  rate.  That  is,  the  swirling  flow  in  annular  fuel  sheet  in  a  rotating  surrounding  gas 
reacts  like  the  inviscid  flow  without  the  surface  tension.  However,  the  effects  of  surface 
tension  are  small  compared  to  the  effect  of  the  swirling  flow  in  annular  sheet.  The  effects 
of  surface  tension  are  significant  in  Reshotko  and  Monnin's  two-wheel  flow.  However, 
the  growth  rate  does  not  continually  increase  as  the  axial  wave  number  increases,  as  in 
Reshotko's  two  wheel  flows  without  surface  tension,  but  begins  to  decrease  regardless  of 
surface  tension.  This  result  is  worthy  of  notice.  Regardless  of  whether  surface  tension 
exists  or  not,  maximum  growth  rate  at  K  =  0.6  is  shown  in  figure  16  and  17.  The 
maximum  growth  rate  occurs  at  small  axial  wave  number  and  high  azimuthal  wave 
number.  Because  of  the  limitation  of  the  computational  capability  on  CRAYC90,  the 
azimuthal  wave  number  beyond  m  =  20  could  not  be  investigated.  It  seems  that  the 
growth  rate  increases  as  m  increases  in  those  figures.  But,  the  difference  of  growth  rate 
between  m  =  0  and  m  =  5  is  larger  than  the  difference  of  growth  rate  between  m  =  5  and 
m  =  1 0.  When  the  extrapolation  beyond  m  =  20  at  K  =  0  is  made  based  on  the  decreasing 
rate  of  the  growth  rate  between  m  =  0  and  m  =  20,  it  is  predicted  from  this  trend  that  the 
growth  rate  at  K  =  0  may  decrease  around  m  =  60  as  m  increases  even  though  the 
computation  of  that  prediction  is  not  available. 
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Notice  the  magnitude  of  growth  rate  when  comparing  the  axial  annular  jet  as 
shown  in  figure  12  and  the  rotating  flow  of  figure  16.  This  tremendous  increase  in 
growth  rate  magnifies  the  impact  of  swirling  flow  of  the  annular  fuel  sheet  on  fuel 
atomization.  In  figure  18,  the  contour  plot  of  growth  rate  is  shown  to  make  figure  16  and 
17  more  understandable.  The  phenomena  that  the  tremendous  maximum  growth  rate 
occurs  at  small  K  such  as  K  =  0.6  shows  that  predominent  instability  due  to  centrifugal 
force  effect  has  an  immediate  impact  on  the  breakup  of  fuel  sheet. 

In  figure  19,  the  effect  of  density  increase  in  annular  fuel  sheet  in  the  rotating 
surrounding  gas  is  investigated.  The  maximum  growth  rate  occurs  at  K  =  0.2.  The 
growth  rate  slowly  decreases  as  K  increases  beyond  K  =  0.2.  This  shows  that  the  effect 
of  density  increase  exist  to  some  extent  over  small  K.  However,  the  effect  of  density 
increase  of  the  annular  sheet  causes  the  shift  of  axial  wave  number  for  the  maximum 
growth  rate  to  take  place  and  the  growth  rate  within  small  K  to  slowly  change.  The 
phenomena  that  the  maximum  growth  rate  occurs  not  at  K  =  0.6  but  at  K  =  0.2  by  density 
increase  of  the  annular  fuel  sheet  tells  that  the  more  severe  breakup  of  the  annular  fuel 
sheet  occurs  at  larger  density  of  the  annular  sheet  than  at  the  typical  density  of  annular 

fuel  sheet.  That  prediction  is  expected  in  terms  of  pitch  angle  tan  6  =  — .  The  inclined 

rk 

degree  of  6  shows  the  intensity  ratio  of  swirling  disturbance  with  respect  to  axial 
disturbance  for  maximum  growth  rate.  The  pitch  angle  is  steeper  when  axial  wave 
number  k  for  maximum  growth  rate  decreases.    One  noticeable  point  is  that  the 
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Magnitude  of  the  maximum  growth  rate  at  K  =  0.6  increases  a  little  as  the  density  of 
annular  sheet  increases  and  growth  rates  increase  large  at  small  K  except  K  =  0.6. 

The  reason  the  maximum  growth  rate  does  not  significantly  increase,  is 
hypothesized  as  a  result  of  the  cohesive  nature  of  increased  density  counteracting  the 
accelerating  effect  of  density  of  annular  fuel  sheet  in  disintegration  by  centrifugal  force. 
The  phenomena  of  the  density  increase  causing  the  growth  rate  to  decrease  has  been 
shown  and  explained  for  the  axial  annular  sheet  jet  in  figure  13.  Such  a  phenomena  is 
thought  to  occur  in  the  similar  characteristic  for  even  swirling  flow  of  annular  fuel  sheet. 
In  figure  20  compared  to  figure  1 8,  the  contour  plot  for  the  effect  of  density  increase  is 
shown  to  help  elucidate  this  concept. 

In  figure  21,  the  effects  of  thickness  decrease  in  annular  fuel  sheet  is  shown.  The 
growth  rate  remarkably  increases  as  the  thickness  decreases.  In  comparison  with  figure 
16  of  typical  thickness,  the  growth  rate  increases  substantially  as  the  thickness  decreases. 
A  similar  phenomena  was  shown  for  the  axial  annular  jet  in  figure  14.  Also,  the 
maximum  growth  rate  occurs  at  K  =  0.2  as  that  occurs  at  K  =  0.2  for  the  density  increase 
of  the  annular  fuel  sheet.  The  pitch  angle  9  for  thickness  decrease  is  steeper  than  6  for 

the  typical  thickness.  It  is  predicted  that  thinner  fuel  sheet  disintegrates  much  faster  into 
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drops.  This  highlights  the  importance  of  presence  of  the  thinner  sheet  for  good 
atomization  quality.  The  annular  fuel  sheet  thickness  may  therefore  be  of  crucial 
importance  to  the  eventual  drop  sizes.  The  contour  plot  of  the  effects  of  thickness  is 
shown  in  figure  22. 
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In  figure  23,  the  effect  of  angular  velocity  increase  in  annular  fuel  sheet  is  shown, 
the  growth  rate  increases  a  little  as  the  angular  velocity  increases  in  comparison  with 
figure  16.  When  angular  velocity  increases  in  annular  fuel  sheet,  the  number  of  azimuthal 
mode(m)  is  increased.  This  phenomenon  is  shown  in  figure  137  of  Van  Dyke  (1980)'s 
album  of  fluid  motion.  The  increase  of  the  number  of  azimuthal  mode(m)  cause  the 
tangential  wave  number  of  figure  5  to  increase.  It  promote  the  breakup  of  the  annular 
sheet  to  many  fine  sizes.  The  reason  the  maximum  growth  rate  does  not  significantly 
increase  is  an  unknown  priori.  The  contour  plot  of  the  effects  of  thickness  is  shown  in 
figure  24. 

For  large  axial  wave  numbers,  the  growth  rates  continually  decrease  as  K 
increases  beyond  K  =  1.  In  figure  25,  the  growth  rate  is  shown  for  the  typical  model  at 
large  K.  The  growth  rate  at  K  =  0  and  m  =  20  or  m  =  15  are  less  than  those  at  K  =  1  and 
m  =  20  or  m  =  15.  On  the  other  hand,  the  growth  rates  at  K  =  0  and  m  =  10  or  m  =  5  are 

larger  than  those  at  K  =  1  and  m  =  10  or  m  =  5.  The  growth  rates  of  m  seem  to  converge 

■ 

as  K  increases.  Therefore,  the  growth  rates  at  m  =  15  or  m  =  20  rapidly  decrease  to 
converge  at  K  =  7  or  larger  K.  In  figure  26,  the  contour  plot  for  the  growth  rates  is  shown 
for  the  typical  model  at  large  K.  In  figure  27,  the  growth  rate  at  K  =  1  for  density  increase 
is  higher  than  the  growth  rate  at  K  =  1  for  the  typical  model  shown  in  figure  25.  Also,  the 
growth  rates  regardless  of  m  converge  at  K  =  9  or  larger.  Therefore,  the  growth  rate 
gradually  decrease  as  K  increases  beyond  K  =  1  and  approaches  9.  The  growth  rate  for 
density  increase  in  figure  27  is  higher  than  the  growth  rate  for  typical  model  in  figure  25 
as  K  approaches  9.  In  figure  28,  the  contour  plot  for  the  growth  rate  is  shown  for  density 
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Figure  24.  The  effect  of  angular  velocity  increase 
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Figure  25.  Typical  model  for  large  K 
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Figure  26.  Typical  model  for  large  K 
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increase.  In  figure  29,  the  growth  rate  at  K  =  0  for  thickness  decrease  is  much  higher  than 
the  growth  rate  at  K  =  1 .  This  is  due  to  the  very  steep  decrease  of  the  growth  rate  as  K 
approaches  1  in  figure  21.  Also,  the  growth  rates  regardless  of  m  converge  at  K  =  2  or 
larger.  Therefore,  the  growth  rate  rapidly  decreases  as  K  approaches  2.  The  growth  rate 
for  thickness  decrease  in  figure  29  is  lower  than  the  growth  rate  of  typical  model  in  figure 
25  as  K  approaches  2.  In  figure  30,  the  contour  plot  for  the  growth  rate  is  shown.  In 
figure  3 1 ,  the  growth  rate  is  shown  for  the  effect  of  angular  velocity  increase  at  large  K. 
The  growth  rates  at  K  =  0  and  K  =  1  are  higher  than  the  growth  rate  of  typical  model  in 
figure  25.  The  growth  rates  of  the  rest  of  m  seem  to  converge  as  K  increases.  In  figure 
26,  the  contour  plot  for  the  growth  rate  is  shown. 

6.4  Helical  Flow  Annular  Jet 

The  investigation  of  the  instability  of  the  annular  fuel  sheet  both  moving  in  the 
axial  direction  and  swirling  in  the  circumferential  direction  in  a  rotating  surrounding  gas 
has  been  done.  Equation  (4.57)  for  m  =  0  and  K  *  0,  eq.  (4.61)  for  m  *  0  and  K  =  0,  eq. 
(4.63)  for  m  *  0  and  K  ^  0  are  solved  for  no  surface  tension.  Equation  (4.58)  for  m  =  0 
and  K  *  0,  eq.  (4.62)  for  m  *  0  and  K  =  0,  eq.  (4.64)  for  m  *  0  and  K  *  0  are  solved  for 
the  effect  of  surface  tension. 

In  figure  33,  the  growth  rate  of  the  typical  model  for  small  K  is  shown.  It  shows 
that  there  seems  to  be  no  difference  between  figure  16  for  the  typical  model  of  swirling 
flow  and  figure  33.  This  tells  that  the  degree  of  instability  of  the  annular  swirling  sheet 
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in  a  rotating  surrounding  gas  is  almost  the  same  as  that  of  the  instability  of  the 
combination  of  annular  axial  jet  and  annular  swirling  sheet  in  a  rotating  surrounding  gas. 
In  figure  34,  the  effects  of  density  increase  are  shown.  In  figure  35,  the  effect  of 
thickness  decrease  is  shown.  In  figure  36,  the  effect  of  angular  velocity  increase  is 
shown.  Comparing  figure  34  with  figure  19  and  figure  35  with  figure  21  and  figure  23 
with  figure  36,  it  is  shown  that  there  is  no  difference  between  the  instability  of  annular 
swirling  fuel  sheet  and  the  instability  of  the  combination  of  annular  fuel  sheet  moving 
axially  and  swirling  circumferentially  in  a  rotating  environment.  In  addition,  even  the 
increase  of  axial  velocity  increase  for  the  combination  of  the  annular  sheet  axially  moving 
and  swirling  circumferentially  does  not  affect  the  growth  rate  as  shown  in  figure  37.  For 
the  annular  axial  jet ,  there  is  increase  of  growth  rate  in  figure  15  when  the  axial  velocity 
in  the  annular  axial  sheet  increases.  However,  this  phenomena  is  not  shown  at  all  in  the 
combination  of  annular  fuel  sheet  axially  moving  and  swirling  circumferentially  in  a 
rotating  surrounding  gas. 

Finally,  it  can  be  said  that  the  swirling  motion  of  annular  fuel  sheet  in  the  rotating 
gas  encroaches  upon  relative  effects  of  annular  sheet  moving  in  the  axial  direction.  It 
causes  strong  centrifugal  instability  to  take  place.  Therefore,  strong  centrifugal  instability 
contributes  to  much  faster  breakup  into  drops  from  the  annular  fuel  sheet  as  soon  as  the 
annular  fuel  sheet  is  exposed  to  the  rotating  surrounding  gas  in  gas  turbine  combustor. 
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Figure  33.  Typical  model  for  small  K 
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Figure  21 .  The  effect  of  thickness  decrease 
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Figure  23.  The  effect  of  angular  velocity  increase 
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Figure  36.  The  effect  of  angular  velocity  increase 
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Figure  33.  Typical  model  for  small  K 
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CHAPTER  7 
CONCLUSION 

Due  to  the  increasing  application  of  airblast  atomizers  in  modern  gas  turbine 
engines,  all  the  variables  which  may  have  some  effect  on  the  process  of  atomization 
should  be  investigated.  So  far  very  little  has  been  done  mathematically  concerning  the 
mechanism  of  disintegration  of  annular  fuel  sheet  exposed  to  the  rotating  surrounding  gas 
in  gas  turbine  combustor. 

The  degree  of  disintegration  of  annular  fuel  sheet  is  due  in  part  to  the  intensity  of 
the  instability  of  the  annular  fuel.  This  process  of  disintegration  is  completed  just  after 
the  nozzle  tip  as  shown  in  the  experiment  photographs.  The  variables  which  effect  the 
instability  growth  rate  are  surface  tension,  density,  fuel  sheet  thickness,  axial  velocity 
and  tangential  velocity. 

7.1  Axial  Flow  Annular  Jet 


1 .  The  instability  growth  rate  decreases  at  a  much  larger  value  of  surface  tension 
than  is  typical  of  liquid  fuel. 

2.  The  growth  rate  decreases  as  the  density  increases. 

3.  The  growth  rate  increases  as  the  fuel  sheet  thickness  decreases. 
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4.  The  growth  rate  increases  as  the  axial  velocity  increases. 

7.2  Rotating  Flow  Annular  Jet 

1 .  The  role  of  surface  tension  is  insignificant.  The  growth  rate  for  rotating  flow 
annular  jet  is  much  higher  than  that  for  axial  flow  annular  jet. 

2.  The  role  of  density  increase  is  to  increase  the  centrifugal  effects  by  the 
increase  of  the  centrifugal  force  of  the  annular  fuel  sheet.  Also,  increased 
density  provides  more  resistance  to  disintegration  of  the  annular  fuel  sheet. 
Therefore,  the  growth  rate  does  not  increase  much  but  does  increase  as  the 
density  increases. 

3.  The  growth  rate  increases  as  the  fuel  sheet  thickness  decreases. 

4.  The  growth  rate  is  slightly  higher  as  the  angular  velocity  increases.  The  reason 
the  growth  rate  does  not  greatly  increase,  is  not  known. 

7.3  Helical  Flow  Annular  Jet 

Items  1,  2,  3,  and  4  in  section  7.2  also  apply.  The  instability  growth  rate  showed 
no  difference  between  rotating  flow  annular  jet  and  helical  flow  annular  jet.  The  axial 
velocity  increase  does  not  change  of  the  growth  rate.  The  growth  rate  for  helical  flow 
annular  jet  is  much  larger  at  small  wave  numbers  than  that  of  the  axial  flow  annular  jet. 
Also,  the  growth  rate  for  the  helical  sheet  increases  with  angular  wave  number  while  the 
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growth  rate  for  the  helical  sheet  increases  only  at  very  small  axial  wave  numbers  and  then 
decreases.  These  results  seem  to  be  enough  to  explain  the  phenomena  of  sudden  fuel 
sheet  breakup  just  after  the  nozzle  tip. 

The  swirling  effect  of  annular  fuel  sheet  in  rotating  environment  has  one  of  the 
greatest  impacts  on  strong  helical  breakup  of  annular  fuel  sheet  in  gas  turbine  combustor. 
The  effect  of  an  annular  fuel  sheet  thickness  decrease  and  angular  velocity  increase  are 
relatively  strong  in  comparison  with  other  effects.  They  are  the  best  candidates  to  be  used 
in  the  design  of  atomizers. 

In  future,  better  asymptotic  expressions  of  Bessel  functions  in  the  dispersion 
equation  are  required  to  obtain  accurately  computational  results  after  solving  the 
dispersion  equation.  Mathematical  formulation  including  nonlinear  higher  order  terms  is 
desirable  to  research  for  a  more  physical  approach  to  actual  breakup  phenomena.  Also, 
instability  analysis  considering  viscosity  is  desirable  to  research  for  the  other  types  of 
atomizers  in  which  viscosity  effect  may  be  important. 


APPENDIX 

COMPUTATIONAL  SOLUTION  METHODOLOGY 


The  dispersion  equations  derived  in  Chapter  4  are  the  two  types  of 
H(iKq)   J(iKq)  p, 

These  equations  are  very  difficult  solve  directly  on  a  supercomputer  since  they  include 
Bessel  functions  and  the  complex  arguments  are  in  Bessel  functions.  It  is  possible  to  get 
the  asymptotic  solutions  for  these  types  of  equations  by  changing  the  dispersion  equation 
into  asymptotic  expression.  When  the  absolute  magnitudes  of  the  arguments  in  Bessel 
functions  are  toward  either  0  or  infinity,  the  asymptotic  expressions  are  available. 
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2^ 


and  Kq^^  =  -Kq+.... 
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When  the  above  asymptotic  forms  are  substituted  into  the  dispersion  equations  in 
Chapter  4,  the  dispersion  equations  are  solvable.  The  disadvantage  of  this  asymptotic 
approach  seems  that  the  asymptotic  form  does  not  depend  on  K  when    Kq  — >  0 . 

Therefore,  the  asymptotic  solution  independent  on  K  is  obtained  when  Kq  — >  0 .  The 

asymptotic  solution  dependent  on  K  is  only  obtained  when  Kq  — >  oo .  For  this  reason, 

the  asymptotic  solution  for  Kq  — >  0  is  regarded  as  the  asymptotic  solution  at  K  -  0. 

Also,  the  asymptotic  solution  for  Kq  — >  oo  is  regarded  as  the  asymptotic  solution  at 

K*  0. 

The  Flowchart  for  the  computer  code  to  validate  the  asymptotic  approach  in 
comparision  with  Reshotko  and  Monnin's  results  is  as  follows: 
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The  Flowchart  for  the  rotating  flow  annular  jet  is  shown  to  investigate  the 
instability  growth  rate  as  follows: 
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